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ABSTRACT 

For  the  flow  of  a  rarefied  gas  of  Infinite  extent 
past  an  Immersed  body  at  large  Knudsen  number,  the  des- 
cription of  the  flow  field  by  the  free-molecule  solution 
Is  only  valid  in  the  Inner  region,  a  region  within  a 
distance  much  larger  than  the  body  dimension  but  much  less 
than  a  mean  free  path  away  from  the  body.   Outside  of  this 
region  the  colllslonal  effects  cannot  be  considered  as  a 
perturbation  even  In  the  lowest  order  approximation.   fb 
secure  a  systematic  method  of  obtaining  high  order  col- 
llslonal corrections,  It  would  be  necessary  to  treat 
adequately  the  outer  region,  the  flow  field  outside  the 
Inner  region.   For  this  purpose  an  asymptotic  theory  for 
obtaining  nearly  free-molecule  solutions  which  are 
uniformly  valid  throughout  the  flow  field  is  proposed  in 
the  present  study  for  the  linearized  Boltzmann  equation 
with  hard-sphere  molecules.   Power-law  molecules  can  be 
treated  in  a  like  manner.   The  theory  is  based  on  a  new 
Integral  formulation  of  the  linearized  Boltzmann  equation. 
An  Inner-and-outer-expan sion  procedure  is  used.   The 
simplification  in  the  inner  region  is  that  the  colllslonal 
effects  are  only  secondary  while  outside  of  the  inner  region 
the  angle  subtended  by  the  body  is  small  so  that  the  flow 
field  is  essentially  that  created  by  a  source  at  the  origin. 
The  form  of  the  source  term  is  determined  by  the  free- 
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molecule  solution  to  the  problem.   The  nature  of  the  flow 

field  in  the  outer  region  varies  from  being  free-molecule 

to  collision-cominated .   The  effects  of  the  lowest  order 

outer  solution  on  the  inner  solutions  are  of  order  l/K 

'    n 

in  two-dimensional  problems  and  of  order  l/K  in  three- 
dimensional  problems. 

For  the  aerodynamic  drag,  it  is  found  that  for  two- 
dimensional  problems 


D  -  D^  M  +  a(^n  -)D.  +  aD. 
I  -M  ■        a '  i      c 


and  for  three-dimensional  problems, 

2    1        2 
D  =  D„  ..   +  aD,  +  (a  in  -)Do  +  a  D^  . 
f  .M      1    ^      a   2      p 

In  two-dimensional  cases,  D-,  can  be  calculated  by- 
integral  iteration,  and  Dp  is  partially  due  to  the  effects 
of  the  lowest  order  outer  solution.   For  three-dimensional 
problems,  D^  is  partially  due  to  the  effects  of  the  lowest 
order  outer  solution  while  D-,  and  Dp  can  be  calculated  by 
the  integral  iteration. 

In  order  to  examine  the  mathematical  nature  of  this 
asymptotic  theory,  a  simple  problem  of  a  two-dimensional 
circular  cylinder  rotating  in  a  rarefied  gas,  which  is  at 
rest  at  infinity,  is  studied.   The  theory  is  applied  to 
this  problem,  using  a  model  equation.   It  is  proved  that 
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the  inner  and  outer  solutions  obtained  by  the  theory  are 
Indeed  asymptotic  to  the  true  solution  throughout  the 
physical  space. 
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I.   Introduction 

I . 1 .   Aerodynamics  of  Rarefied  Gas  Flows 

For  a  gas  flow,  when  the  mean  free  path  A  of  the  gas  Is 
not  negligibly  small  In  comparison  with  the  characteristic 
length  of  the  problem,  the  Navler-Stokes  stress  vs.  rate-of- 
straln  relations  and  the  Fourier  heat  flux  law  no  longer  are 
valid.   Consequently,  the  Navler-Stokes  equations  of  gas 
dynamics  have  to  be  abandoned  In  favor  of  a  kinetic  theory 
treatment.   The  experimental  results  on  the  low  pressure  gas 
flows  through  tubes  obtained  by  Knudsen  (1909)  and  Gaede  (1913) 
clearly  demonstrated  this  failure  of  the  Navler-Stokes  equa- 
tions.  The  length  L  of  the  problem  may  be  taken  to  be,  for 
example,  the  characteristic  dimension  of  the  body  present 
In  the  flow,  the  boundary  layer  thickness,  or  the  thickness 
of  a  shock  wave;  Its  choice  depends  on  the  aspects  of  the 
problem  that  is  of  Interest.   A  kinetic  theory  treatment  of 
gas  flow  problems  is  not  available  in  full  generality  because 
of  the  severe  mathematical  difficulties  involved  in  constructing 
solutions  to  the  Boltzmann  equation,  the  basic  equation  of  the 
kinetic  theory  of  neutral  gases.   However,  for  a  gas  flow  of 
infinite  extent  past  an  Immersed  body,  when  the  characteristic 
dimension  L  of  the  body  is  very  small  in  comparison  with  the 
mean  free  path  A  of  the  gas,  i.e.,  when  a  =  L/A  «  1,  one 
method  of  constructing  an  approximate  solution  to  the  flow 
problem  is  to  neglect  the  Interaction  between  the  molecules 
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re-emltted  from  the  body  and  the  incident  stream,  cf .  Epstein 
(1924)  and  Tslen  (19^6).   The  solution  so  obtained  is  called 
the  free-molecule  solution.   The  simplifying  feature  of  such 
flow  problems  is  that  most  of  the  molecules  re-emitted  from 
the  body  do  not  collide  with  the  molecules  in  the  incident 
stream  until  very  far  away  from  the  body.   This  justifies  the 
separate  consideration  for  the  re-emitted  molecules  and  the 
molecules  in  the  incident  stream  in  the  lowest  order  approxima- 
tion.  This  type  of  problem  is  of  practical  interest  in  con- 
nection with  flights  at  very  high  altitudes  where  the  mean 
free  path  of  the  atmosphere  is  very  large.   The  free -molecule 
solutions  have  been  found  in  good  agreement  with  experiments, 
cf.  Stalder,  et  al .  (1950). 

1.2 .   Colllsional  Effects  and  Nonuniform  Validity  of 
Free -molecule  Solutions 

A  number  of  Improvements  over  the  free-molecule  solution 
have  been  suggested  in  order  to  take  into  account  the  colllsional 
effects,  e.g.  by  Lunc  and  Lubonskl  (1956),  Baker  and  Charwat 
(1958),  and  Willis  (I958).   Willis  has  transformed  the 
Boltzmann  equation  into  an  integral  equation  somewhat  similar 
to  the  Voterra  type,  and  the  method  of  solution  he  has  proposed 
is  an  iteration  procedure  for  this  Integral  equation  using  the 
free -molecule  solution  as  the  zeroth  iterate.   Taking  the 
Boltzmann  equation  as  the  basic  framework,  Willis  has  been 
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able  to  eliminate  the  unnecessary  assumptions  and  restrictions 
In  the  previous  works.   Solutions  that  take  Into  account  the 
colllslonal  effects  for  L/A  «  1  are  called  nearly  free- 
molecule  solutions.   In  essence,  Willis'  method  reduces  the 
procedure  of  finding  nearly  free-molecule  solutions  for  linear 
as  well  as  non-linear  rarefied  gas  flows  to  quadrature  compu- 
tations.  It  is  claimed  that  more  accurate  correction  for 
colllslonal  effect  can  be  obtained  by  going  to  higher  itera- 
tions . 

However,  it  should  be  noted  that  the  description  of  the 
flow  field  by  the  free-molecule  solution  is  only  valid  in  a 
region  within  a  distance  much  larger  than  the  body  dimension 
L  but  much  less  than  the  mean  free  path  A  from  the  body.   This 
region  may  be  called  the  inner  region.   Outside  of  this  region 
the  colllslonal  effects  cannot  be  considered  as  a  perturbation 
even  in  the  lowest  order  approximation,  and  as  a  result  the 
free-molecule  solution  ceases  to  be  a  proper  description  of 
the  flow  field.   Because  of  this  nonuniform  validity  of  the 
free-molecule  solution,  any  iteration  procedure  based  on  the 
free -molecule  solution  as  the  zeroth  iterate  will,  after  a 
certain  Iterative  state,  have  to  fail  to  give  correct  account 
for  the  colllslonal  effects. 


-  5 


1.3-   Purpose  of  the  Present  Study 

1 

It  Is  thus  clear  that  a  systematic  method  of  obtaining 
high  order  corrections  for  the  colllslonal  effects  cannot  be 
achieved  based  on  the  free-molecule  solution  only,  and  that 
to  arrive  at  such  a  method  It  Is  necessary  to  treat  adequately 
the  outside  region.  I.e.  the  part  of  the  physical  space  out- 
side of  the  Inner  region.   In  this  connection,  a  question  of 
practical  Importance  arises  as  to  whether  the  first  and 
second  Iteration  results  for  aerodynamic  forces  obtained  by 
Willis  '  method  as  a  higher  order  approximation  will  be 
affected  by  the  nonuniform  validity  of  the  free-molecule 
solution. 

It  Is  the  Intention  of  the  present  study  to  give  an 
asymptotic  theory  for  obtaining  nearly  free-molecule  solutions 
which  are  uniformly  valid  throughout  the  physical  space  and  to 
answer  the  question  mentioned  above.   This  asymptotic  theory 
Is  based  on  a  new  Integral  equation  formulation  for  the 
linearized  Boltzmann  equation.   The  Integral  equation  Is 
Predholm  In  c_,  the  molecular  velocity  and  Is  derived  In 
Chapter  II.   The  theory  essentially  uses  an  Inner  and  outer 
expansion  procedure  and  Is  carried  out  for  hard -sphere 
molecules  In  Chapter  III.   Power-law  molecules  with  a  finite 
cut-off  can  be  treated  In  a  like  manner.   In  the  Inner  region 
the  simplification  is  that  the  colllslonal  effects  are  only 
secondary,  while  in  the  outer  region  the  angle  subtended  by 
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the  body  Is  very  small.   It  Is  found  that  the  outer  expan- 
sions are  governed  by  a  sequence  of  integral  equations  each 
of  which  Is  Independent  of  the  parameter  a.   In  particular, 
the  lowest  order  outer  expansion  corresponds  to  a  flow  field 
created  by  a  source  at  the  origin,  the  form  of  this  source 
term  being  determined  by  the  free -molecule  solution  to  the 
problem.   For  2-dlmenslonal  problems  the  first  order  Inner 
solution  Is  not  affected  by  the  outer  flow  field  and  there- 
fore can  be  computed  by  Willis'  integral  iteration  procedure. 
In  ^-dimensional  problems  this  is  true  for  the  first  and  the 
second  order  inner  solutions.   The  effects  of  the  outer  flow 
field  are  present  in  the  higher  order  inner  solutions. 
The  expressions  for  the  drag  are  found  to  be,  for 
2-dimensional  flows 


D  =  D^^j^_  +  a  -«n(^)D^  +  aDg 


and,  for  3-dimensional  flows 

D  =  D^_^_  +  ccD-^  +   a^in(^)D2  +  ^^^03  . 

In  2-dlmenslonal  cases,  D-,  can  be  computed  by  Willis' 
integral  iteration  method,  and  this  is  true  for  D-,  and  D^ 
in  3-dimenslonal  cases.   To  calculate  the  next  order  collisional 
corrections,  the  lowest  order  outer  solution  has  to  be  calculated 
numerically. 
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In  order  to  Investigate  the  mathematical  correctness 
of  this  asymptotic  theory,  the  simple  problem  of  a  circular 
cylinder  rotating  In  an  Infinite  gas  which  Is  at  rest  at 
Infinity  Is  worked  out  In  chapters  rv  and  V  as  an  Illus- 
trative example  In  considerable  detail  based  on  the  model 
equation  proposed  by  Bhatnagar  et  al .  (1954)  and  Welander 
(1954).   For  this  problem  It  can  be  proved  with  mathematical 
rigor  that  the  nearly  free -molecule  solutions  obtained  by 
the  theory  are  Indeed  uniformly  valid  approximate  solutions 
to  the  problem. 
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II.  Boundary  Value  Problems  for  the  Linearized  Boltzmann 
Equation 

I I . 1 .  Boltzmann  Equation  In  Orthogonal  Curvilinear 
Coordinates 

In  an  orthogonal  curvilinear  coordinate  system  (x.^x^^x^) 
with  its  metric  tensor  g. .  =  g.h-  .,    the  dynamical  state  of  a 
classical  particle  can  be  identified  with  6  numbers  (x-,,Xp,x^; 
X-,  ,x  ,x^)  where  the  dot  denotes  total  differentiation  with 
respect  to  time  t.   The  summation  convention  Is  not  adopted 
in  the  present  work.   An  alternative  description  of  this 
dynamical  state  is  (x-,,Xp,x^;  |  ,^„,£,  ),  or,  in  a  short-hand 
form,  (x.,^.)  where  ^.  =  g.x..   The  state  at  time  t  of  a  gas 
consisting  of  monatomic  molecules,  considered  as  classical 
particles,  is  completely  specified  by  the  number  of  molecules 
having  their  dynamical  states  in  the  range  between  (x. ;^ . ) 
and  (x  +  Ax. ;  I.  +  A^ . )  at  time  t.   The  number  may  be  denoted 
by  P(x. ;| . ;t )Ax. A^ .  where  P(x.;^.;t)  is  called  the  distribu- 
tion function.   The  change  of  P  along  a  molecular  trajectory 
is  due  to  the  collisions  in  the  gas.   This  statement  may  be 
written  as 


dP  _  5P 
dt    6t 


)      ,  (2.1) 


coll 


where  the  R.H.S.  term  stands  symbolically  for  the  colllsional 
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effects  on  F,  and  may  be  evaluated  In  terms  of  F  If  the 
assymptlons  of 

(1)  binary  collisions, 

(2)  molecular  chaos, 

(3)  finite  molecular  cross  section,  and 

(^)  slowly  varying  F  over  a  molecular  dimension 
are  made,  which  will  always  be  assumed  to  be  the  case  through- 
out this  work.   Under  these  assumptions,  equation  (2.1)  is 
the  renowned  Boltzmann  equation. 

The  L.H.S.  of  equation  (2.1)  stands  for  the  limit 

F(x^  +  5x^;  i^+^i^;    t+6t)  -  F(x^;^^;t) 

lim ■ ■ 

6t  — >  o  5t 

taking  along  the  molecular  trajectory,  i.e., 
dF   SF  _^  ^  ^i  SF   ,  ^  A      aP 

The  factor  ^.  must  be  evaluated  in  terms  of  the  state 
variables  (x. ;  i.)-  This  may  be  done  most  conveniently  by 
using  the  Lagrange  equations 


d  /SL  \    bL     _  ^ 
dt  ^M7^  ~  ^S^  -  ^    ' 


where  L  =  T-V,  T  Is  the  kinetic  energy  of  a  particle  and  V 
is  its  potential  energy. 
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5 


L  =  IT   ?J  ?f.  I  -  V 


k=l 


2  ^k 


It  then  follows  Immediately  from  the  Lagrange  equations 


that 


d^ 


1    1 


dt     g. 


1 

k= 


3 


(^ 


^Si,      Sg. 


K   /S      ~K     £      -"In      ^1   BV 


^ 


— T-  g,  ^  k  dx.     1  "Sx 
=1   ^k       1 


k 


m 


^5r; 


On  substituting  this  result,  the  Boltzmann  equation 
becomes 


^  4,       ' 


"St"     ^^— ^  g.  dx.  4—rr    g. 

4,=1  '^'l    1    1  =  1  ^1 


3 


^?(? 


Sgv       ^g 


=k 


k=l  ^k 


'k  dx_.    i  "5x 


1 )  _  fl  Sv 

,  '       m  dx. 
k        D 


dF 


5F 
6t 


coll 


(2.2) 


This  is  the  non-linear  Boltzmann  equation  in  orthogonal 
curvilinear  coordinates.   In  the  case  of  cylindrical  coordi- 
nates, (x^^x^^x^)  =(r,0,z)  and  (g^,g2.gj)  =  (l,r,l).   The 
Boltzmann  equation  may  be  written,  in  the  absence  of  force 
field,  explicitly  as 


SF  ^t       BF  ^^0  SF  ,  £   SF  ,  ^0  dF    ^e^r  BF 
^^'^r  ^+~  M+^z  ^  +  ' 


^ 


^ 


5F 

6t 


coll 


(2.3) 


For  gases  of  slight  nonunif ormlty  about  an  equilibrium 
state  of  number  density  n  ,  temperature  T  and  uniform 
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velocity  ^^_o-   uy2RT   ,  It  Is  expedient  to  Introduce 

c.  =  I  ./)/2RT  ,    for  1  =  v,9,z,    and  the  function  f(c  ,Cq,c  ; 
1    1  '    o  '  r  y   z 

r,9,z;    t)  as  a  dependent  variable  through  the  relation 


P  =  Fj^(l  +  f) 

1    ^/2    r        2) 
where   P^  =  n  (g     )    expj-(c_  -  U)  [  •   In  particular,  for 

a  uniform  flow  past  a  body,  n  ,  T   and  u^wlll  be  taken  as 
the  number  density,  temperature  and  velocity  of  the  undis- 
turbed flow. 

On  neglecting  the  terms  quadratic  In  f,  equation  (2.3) 
becomes.  In  2-dlmenslonal  steady-state  problems, 

2 
^   Sf  _^  ^0  Sf  _^  °0  hf  ^0^r  bf  fT        „„> 

C    ^^  +  -rpr  +  •'s -  ^r =  a((P   -  Vf  ) 

r  dr    r  d0    r  dc      r   dCn 


'0 


where,  for  hard-sphere  molecules, 


2t 
a  =  n  a   L 
o 


i(c;r,0,z)    =  ///  K(c,c-^)f  d^c^ 


K(c_,c_-j^)    =   K-|^(£,£-^)    +   K2(c_,c_-j^) 

r-  ,         (       ^  2.   I  (2-5) 

K^(c_,c_^)    -  - /tT  I  £-£i  I  exp    -    (c,)      > 


/—       1  11  2(.C-C,   j  C|«  c 

K2(c,c^)    ^2IT[- exp    -  j^   |c-c_^' 


I  £"£i  I  (      "  "^  •2'  I  £-£]_  I 

.-Ir  C 


v(c)   =    ^  (cj   [  S  expC-c"-)  -,(f  ^pC-^^JJx)(2ci  l)] 
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where  L  is  the  characteristic  length  of  the  body  used  to 
non-dlmenslonallze  equation  (2.7),  and  a  is  the  reciprocal 
of  the  Knudsen  number.   For  the  derivation  of  (2.5)  and  also 
the  corresponding  results  for  power-law  molecules,  cf.,  e.g., 
Grad  (1963)  or  Waldmann  (1958).    The  molecules  described  by 
the  function  f  will  be  called  the  n on -Maxwell Ian  molecules 
and  those  described  by  the  uniform  Maxwelllan  F^  will  be 
called  Maxwelllan  molecules.*   A  molecule  Is  said  to  be  a 
c-molecule  If  Its  velocity  Is  In  the  range  (c_,  dc_) ;  otherwise 
It  Is  said  to  be  a  non-c -molecule .   v(g_)  Is  the  number  of  col- 
lisions encountered  by  a  c-molecule  during  the  time  Interval 
L//2RT  .   The  local  production  ^(c;r,0)  Is  the  net  effect  on 
the  number  of  non -Maxwell Ian  c-molecules  due  to  (l)  collisions 
encountered  at  (r,0)  by  non -Maxwell Ian  non-c -molecules  and 
(11)  collisions  encountered  at  (r,0)  by  Maxwelllan  c-molecules. 
It  should  be  noted  that  the  local  production  (^  may  take  nega- 
tive values  because  the  effect  of  collisions  of  the  second 
type  Is  always  negative  and  may  surpass  that  of  the  first  type 
in  magnitude.   The  effects  of  the  collisions  of  the  first  and 
the  second  type  are  represented  by  the  actions  of  the  kernels 
K-]  and  Kp  respectively. 


■X- 

This  should  not  be  confused  with  the  usual  reference  of 
Maxwelllan  molecules  to  the  inverse  fifth  power  law  mole- 
cules. 
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II. 2.   Boundary  Conditions. 

Two  types  of  boundary  conditions  arise  In  the  problem 
of  a  uniform  flow  past  an  Immersed  body.   They  are  the 
boundary  condition  at  Infinity  and  that  at  the  body  surface 

The  boundary  condition  at  Infinity  Is  slmplfled  by  the 
use  of  the  function  f(c;r,0)  as  the  dependent  variable  and 
can  be  expressed  simply  as 


11m  f(c;r,e)  =  0,    (c^  <  O) .  (2.6) 

r  ->  00 


That  Is,  the  Incoming  distribution  at  oo  is  Pj^,. 

Boundary  conditions  of  the  other  type  are  more  complicated. 
Their  formulation  requires  knowledge  of  the  Interaction  between 
the  gas  molecules  and  the  solid  surfaces  and  of  the  condition 
of  the  surface. 

If  D(r,e)  Is  an  abstract  function  which  characterizes 
all  the  solid  properties  of  the  body  surface  at  (r,e),  this 
boundary  condition  may  be  put  In  the  following  form,  using  the 
scattering  probability  distribution  P(c_,c_-j^;  D(r,0)), 

F(c;r,0)  =  fff  P(c,c^;  D(r ,6) )P(c^;r ,0)d^c^   (2.?) 


Cni<  0 


for   c^^  >  0, 


where  c   is  the  component  of  the  molecular  velocity  in  the 
direction  of  the  outward-drawn  normal  n  to  the  body  surface 
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at  (r,0).   And,  in  terms  of  the  perturbation  distribution 
function  f,  equation  (2.7)  becomes 


f  (c_;r,e)  =  -2c. U  +  e^ 


J  J 


2 


e   ^(2c.  -U)?  d^c.        (2.8) 


-1  - 


c,  <  0 


'n 


2 
+  e^ 


1 


e   ^  f(c^;r,0)P  d^c^  , 


c   <  0 
"1 


for   c  >  0  . 
n 


The  surface  condition  D(r,e)  may  conceivably  be  determined 
by  the  solid  structure,  the  surface  temperature  and  the  degree 
of  adsorption  of  the  gas  on  the  surface.   The  surface  tempera- 
ture is  dictated  by  the  conservation  principles.   For  example, 
if  the  surface  is  locally  adiabatic,  the  surface  temperature 
must  maintain  a  balance  of  energy  flux  carried  by  the  incoming 
molecules  and  the  re-emitted  ones  at  the  surface.   The  adsorp- 
tion of  gas  on  the  surface  will,  in  general,  depend  on  the 
pressure  exerted  by  the  gas  at  the  surface,  if  no  special 
effort  is  made  to  control  the  adsorption.   A  theoretical  con- 
sideration of  all  these  phenomena  Is  a  very  involved  matter. 
Indeed,  even  if  the  surface  is  a  clean  crystal  surface  at  a 
known  temperature,  a  theoretical  treatment  will  have  to  compute 
the  quantum  mechanical  interaction  between  the  incoming  molecules 
and  the  lattice  surface,  which  is  Itself  a  many-body  problem. 
Thus,  it  would  seem  more  reliable,  at  the  present  time,  to  obtain 
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the  scattering  probability  distribution  P(c_,c_-,;  D(r,0)) 
from  molec-ular  beam  experiments.   For  a  wide  class  of 
engineering  surfaces,  e.g.,  polished  metal  surfaces,  it  has 
been  found  that  diffuse  re-emission  is  a  good  approximation, 
cf.  Estermann  (i960).   For  diffuse  re-emisslon,  the  distribu- 
tion function  of  the  re-emitted  molecules  is  absolute  Max- 
wellian  at  the  temperature  of  the  surface.   The  treatment 
will  further  be  simplified  by  restricting  it  to  locally 
adiabatic  solid  surfaces,  in  which  case  the  constants  in  the 
re-emission  distribution  function  are  uniquely  determined  by 
the  principles  of  conservation  of  energy  and  matter.   As  a 
result,  the  scattering  probability  distribution  P(c_,c_-,;  D(r,0)) 


is  found  to  be 


c^c^       2 
I  (5  -  ='  -  =1  +  -T^K^^~'  (2-9) 


and  the  boundary  condition  becomes 


f(c;r,0)  =  -2c_.U  -  ^  (l  +  ^)U  cos  6^  (2.10) 


TT 


.2  2  2 


e   ^  c^  (^-c^-c^  +-^)f  d^c-^ 


c^^<  0 

for  c   >  0,  , 

n 

where  9^   is  the  angle  between  U  and  n.   For  the  derivation  of     ] 


equation  (2.9),  see  Appendix  I 
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II. 3-   Boundary  Value  Problems  and  Integral  Equations 

A  new  integral  equation  formulation  of  the  boundary 
value  problems  for  the  linearized  Boltzmann  equation  will 
be  developed  in  this  section.   This  formulation  is  not 
limited  to  any  particular  values  of  the  KJnudsen  number; 
however,  it  is  of  particular  convenience  for  an  asymptotic 
theory  to  be  proposed  later  for  nearly  free-molecule  flows. 
Only  two-dinensional  problems  are  considered  here,  while  the 
treatment  of  three-dimensional  problems  can  be  carried  out 
without  conceptual  difficulties  in  a  similar  fashion. 

Consider  a  uniform  flow,  with  undisturbed  number  density 
n  ,  temperature  T  and  velocity  U  )/2RT   ,  past  a  two- 
dimensional  cylinder  whose  cross-sectional  contour  is  convex 
and  defined  by  the  equation 

r  =  (?(0)  (2.11) 

The  coordinate  system  (r,0)  is  so  taken  that  the  piece- 
wise  differentiable  function  Ik^{Q)    is  single-valued.   The 
length  scale  in  the  problem  is  such  that  the  maximum  of  R(0) 
Is  equal  to  unity,  i.e.. 

Max    (p(0)  =  1  . 
0^  0^2Tr 

When  It  is  necessary  to  emphasize  the  components  of  c. 
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f(c  ^c„,c  ;  r,G)    or  f(c  ,3,0;  r,G)   will  be  used  Instead  of 
f(c;r,0)  where  (c  ,s,0)  are  the  cylindrical  coordinates  of  c_. 
The  mathematical  problem  for  this  two-dimensional  flow  is 
then  to  find  a  function  f(c;r,0)  satisfying  the  following 
set  of  equations: 

c,  ^  +  ^  +  ^  1^  -  ^^  1^  =  -(^  -  vf)       (2.12a) 
1  dr    r    r   dc      r   dc^     \:ii     /        \  j 

r  0 

f{^;r,9)    =  ///  k(c,^^)f(c^;r,0)d^c^  (2.12b) 


f  (c,r,0)  — >  0   as    r  — >  «>  ,   for   c  <  0       (2.12c) 

f(c;r,0)  =  -2c.U  +  -^  (l  +  ^)U  cos  0^  (2.12d) 

-c2  2  2 

e   ^  c^  (5  _c2+c^+-^)f(c^;r,0)d^c^ 


IT 


J 


c^  <  0 


'n 


1 


for  c^  >  0,  r  =(^(0) 


In  obtaining  equation  (2.12d),  the  diffuse  boundary 
condition  of  locally  adlabatic  body  surface  has  been  used. 
Boundary  conditions  of  more  general  type  can  be  treated  in  a 
like  manner  for  linearized  problems. 

The  characteristic  equations  for  the  differential 
operator  of  equation  (2.12a)  are 

dr:_  r'd0'    ^'^^r  _  ^'^^0 

.       •    ~    t2   "   '  ,   ' 
r      Cg      Cq  c^Cq 
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dc;  =  0  , 


or,  in  their  Integrated  forms 

c  =   const   =  c  (2.13a) 


'2  '2  2 

c^  +  c    =  const  =  s  (2.13b) 


r 


I  1 


r  Cq  =  const  =  rs  sin  <!>  (2.13c) 


+  *  -  sin  ^(^  sin  *)  ■  (2.13d) 


For  given  constants  c  ,  s,  <t>,  r,    and  6   equations  (2.13) 
define  a  curve  In  the  5-space  (c_',r',e').   This  curve.  In  fact. 
Is  a  molecular  trajectory  passing  through  the  point  (r,0)  In 
the  physical  (r,e) -plane  at  an  Inclination  0-^  with  the 
x-axls  (see  Fig.  l).   The  z-component  of  the  molecular  velocity 
Is  c   and  the  speed  In  the  (r, (9) -plane  Is  s.   It  may  be  noted 
that  equation  (2.13d)  applies  only  to  an  Interval  of  a 
molecular  trajectory  along  which  r'  varies  monotonlcally . 
Along  such  a  part  of  the  molecular  trajectory  equation  (2.12a) 
may  be  Integrated  In  terms  of  ^,  using  equations  (2.13)  with  r' 
as  the  Integration  variable.   The  resultant  equation  Is 


-^('^.sAj^B)  -  {(^'.^'0') 


^   Hfo  c;^-t; 


r-<=L  (2.14) 

/  $(£>^  /^^^ 
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in  which  c_'  and  0'  are  functions  of  x,  <i> ,    6,    and  r'  as 
defined  by  (2.13) • 

Equation  (2.l4)  can  be  made  more  explicit  If  a   Is 
chosen.   For  this  purpose,  the  velocity  space  (c  ,s,0)  at 
(r,0)  Is  divided  Into  three  regions.   These  regions  and  the 
choice  of  a   In  each  of  them  are : 

(l)  Region  I,  in  which  ir/2   <  |  *  |  <  tt.   The  choice  of  a 
in  this  region  is 


a  =  00. 


(ll)  Region  II,  in  which  the  velocity  vector,  if  reversed 
in  direction.  Intersects  the  body  at  some  point  (R(0u^),0.  )• 
If  -©-,  and  (h)   define  the  boundary  rays  of  this  region,  it 
may  be  characterized  by  -©-,  <  <i>   <®  p-   The  choice  of  a   is 

a  =  R^, 

where  R,  =  R(0,  ),  and  R,  is  also  a  function  of  4> ,  r,  and  0. 

(ill)  Region  III,  the  rest  of  the  velocity  space,  in  which 
-7r/2  <  <t>  <  -B-|_  or  (h).  <  *  <  Tr/2  •  The  choice  of  a  in  this 
region  is 


a  =  r  sm  <p  . 

The  case  when  ©>.  or  0^    is  greater  than  ir/2   can  be  handled 
without  any  essential  change  in  the  method. 

With  these  choices  of  a,  equation  (2.14)  becomes,  in 
Region  I 
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In  Region  II, 


(2.15b) 


and  In  Region  III, 


.   o^f    /(->''"       "-^   (2.15c) 


In  this  region  r'  runs  through  a  minimum  r  sin  4). 

For  the  geometry  of  the  three  regions,  see  Fig.  2.   The 
physical  meaning  of  the  terms  In  equations  (2.15)  Is  quite 
clear.   The  non-Maxwelllan  c-molecules  at  (r,e)  can  be  clas- 
sified Into  two  groups,  namely  (l)  caused  by  collisions  at 
(r,0)  of  molecules  coming  directly  from  the  body  surface  and 
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(ll)  by  collisions  at  {r,6)    of  molecules  coining  from  a  point 
(r',0'),  on  the  molecular  trajectory,  where  they  suffered 
their  last  collision.   In  equations  (2.15),  the  terms  involving 
^  represent  the  second  group  of  molecules  and  the  known  func- 
tions in  (2.15b)  represent  the  first  group.   The  exponential 
factors  are  the  probabilities  of  a  c-molecule  to  travel 
without  suffering  any  collisions  from  (r',0')  to  (r,0)  and 
the  denominator  in  the  integrals  is  a  factor  which  relates 
the  differential  length  dr '  to  the  corresponding  differential 
length  along  the  molecular  trajectory. 

The  boundary  condition  (2.12c)  for  f  at  r  =  «>  can  be 
translated  into  a  similar  condition  for  ^,  using  expression 
(2.15a).   For,  according  to  the  mean-value  theorem  in 
integral  calculus,  expression  (2.15a)  may  be  written  as 

where  ^^^^^  is  the  value  of  \{'z  ^,z  ,<!> '  ;    r',^')  taken  at  s 


ome 

^m^or^^  ^  <  ^r^^ov,  <  ~-   ^^  then  follows  that  (2.12c)  is  true 
if  and  only  if 

11m  f  (c^;r,0)  =  0  ,    for  c  <  0  . 
r— ^00 


In  all  the  regions  I,  II,  and  III,  f(c_;r,0)  can  be 
expressed  in  terms  of  known  functions  and  Integrals  Involving 
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^(c_;r,0).   Substituting  these  expressions  Into  (2.12b),  which 
defines  the  function  ^,  gives  rise  to  a  linear  Integral  equa- 
tion for  ^.   If  the  dependence  of  ^  on  the  parameter  a  Is 
expressed  explicitly,  this  Integral  equation  may  be  written 
In  the  following  short -hand  form: 


^(c_;r,0;a)    =  ^  (c_;r,0;a 


(2.16) 


+  aB-^[ 


r  ( 

^  T,dr']    +   aB2[      ^  T,  dr  '  ] 


r  r 

r 


+  al^[  1    $_  T^dr'  ]    +  a1^ 

R^, 


T^dr' ] 


r  slnl4',  I 


+  al^[ 


\  T^dr'  ]    +   a\^  \  "Y^^V  ] 


r  slnlcf'ii 


where  \  must    satisfy  the   condition 


11m     ^(c_;r,e;a)    =   0        for        c^  <    0 
r  — >oo 


(2.17) 


The  full  expressions  of  the  function  ^  and  the  opera- 
tors B's  and  I ' s  In  the  equation  (2./6)  are  as  follows: 

^(c;r,e;a)  =I,[5,{-2C,U  tf  0.c,^jUa>se4  expf-f  (^cos^, '/^^?^)]]  ^ 


«3  ®j. 


O  -(© 
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b^Ie] 


-00     0    -&, 


^     ^     ®z 


l2[g] 


l3[g] 


;        'J        L  -JY      ir/z 


-bO  0 


Tt    = 


l1<^AS 


^.'11 


^2    = 


y^r=^  .»rf-  ^  f/^^T^Iv;^  +  ri^.^;;  i 


T^,    = 


where   R,     Is   the    sol-utlon   to  the    simultaneous   equations 


\  =   K(e,)  , 


e^   =    0  +    <t)^    -    sin   '^  [^  sin   <1>-^) 


and  in  the  Integrands,  the  arguments  of  the  function 
^(c   ,s.,<t)';  r',0';a),   1  =  0,1,  are  understood  to  be 
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-1  V 

4) '  =  sin   ( — r  sin  0  .  ) 


+  0.  -  4)  ' 

1 


Despite  the  complicated  appearance  of  equation  (2.l6), 
the  physical  Interpretation  Is  rather  apparent.  On  the  R.H.S. 
of  (2.l6),  the  first  three  terms  represent  the  contribution 
to  ^  due  to  molecules  coming  directly  from  the  body  surface 
and  the  rest  of  the  terms  that  due  to  molecules  coming  from 
their  last  places  of  collision  in  the  gas. 

The  system  of  equations  (2.12)  has  now  been  case  into 
the  system  of  equations  (2.l6)  and  (2.17).  In  other  words, 
the  mathematical  problem  of  all  the  boundary  value  problems 
for  the  linearized  Boltzmann  equation  is  to  find  a  function 
^(c_;r,0;a)  which  satisfies  equation  (2.l6)  under  the  condi- 
tion (2.17). 

It  is  worthwhile  to  point  out  certain  features  of  this 
integral -equation  formulation.   In  the  system  of  equations 
(2.l6)  and  (2.17)  there  is  only  one  single  Integral  equation 
while  in  the  system  (2.12)  there  are  integro-diff erential 
equation  and  more  complicated  boundary  conditions.   Furthermore, 
equation  (2.l6)  is  Predholm  in  c_;  i.e.,  its  domain  of  integration 
is  independent  of  the  molecular  velocity  c_.   It  is  perhaps  advan- 
tageous in  certain  problems  to  use  $  as  the  dependent  variable 
since  ^  is  a  continuous  function.   It  the  next  chapter  an  asymp- 
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totlc  theory  for  nearly  free -molecule  problems  can  be 
developed,  based  on  equations  (2.l6)  and  (2.I7). 
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Ill-  A  Uniformly  Valid  Asymptotic  Theory 

III.l.  Approximate  Solutions  and  Asymptotic  Expansions 

A  number  of  concepts  and  terms  are  Introduced  here  to 
facilitate  the  asymptotic  analysis  to  be  carried  out  In  the 
later  sections.   First,  consider  the  equation 

L(x,  u(x,a);  a)  =  0   for   x  £  D  ,  (3-1) 

which  Is  assumed  to  define  a  well-posed  problem  with  a 
parameter  a.   By  this  Is  meant  that  there  exists  a  unique 
solution  u(x;  a)  to  {^.l),    uniformly  bounded  In.D.   A  function 
v(x;  a)  Is  said  to  be 

(l)  an  approximate  solution  to  equation  (3-1)  within  an 
absolute  error  e  If 

|u(x;  a)  -  v(x;  a)|  <  e     for  all   x  6  D       (3-2a) 
or 

(11)  an  approximate  solution  to  equation  (3-1)  within  a 
relative  error  e.  If 

|u(x;  a)  -  v(x;  a)|  <  £|u(x;  a)|   for  all   x  6  D.  (3-2b) 

For  a  £  (0,a  )  and  a  «  1,  approximate  solutions  may  be 
constructed  from  knowledge  of  the  nature  of  u(x;  a)  In  the 
limit  a  — >0.   For  example,  suppose  two  sequences  of  functions 
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a  (a)  and  u  (x),  n  =  0,  1,...,  N,  can  be  found  such  that  for 
k  =  0,  1,  . . . ,  N-1, 

lim   ^+1  ^   -  0  (3.3) 

a— >  0    k^  ^ 

k 
u(x;a)  -  ^^  a^(a)u^(x) 

llm  2z^- =  u,^,(x)  (J.4) 

a->0        ^k+l^°^^ 

k 

The  functions  J/^     o    (a)u  (x)  are  said  to  be 

n=o 

(l)  a  k-th  order  asymptotic  expansion  for  u(x;  a)  If 
u,  -,  (x)  Is  uniformly  bounded  In  D,  or 

(11)  a  k-th  order  uniformly  asymptotic  expansion  for 
u(x;  a)  If  u,  .  (x)/u  (x)  Is  uniformly  bounded  In  D  by  M,  ,. 

These  may  be  designated  as  u(x;  a)  ^^  Z   o    (a)u  (x)  for 
an  asymptotic  expansion  and  u(x;  a)  ^^  Z   o    (a)u  (x)  for  a 
uniformly  asymptotic  expansion.   It  can  be  shown  that  Indeed 
an  asymptotic  expansion  for  u(x;  a)  Is  an  approximate  solution 
to  equation  (3-1)  within  some  absolute  error  e  «  1,  and  a 
uniformly  asymptotic  expansion  for  u(x;  a)  Is  an  approximate 
solution  to  (3-1)  within  some  relative  error  e  «  1.   To  do 
this  for  the  asymptotic  expansions  Is  trivial.   In  what 
follows  this  will  be  verified  for  uniformly  asymptotic  expan- 
sions. 

This  Is  easily  seen  for  k  =  0;  In  fact,  given  0  <  e  <  i, 
it  is  always  possible  to  find  a   such  that  for  a  6  (O,  a    ) 
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u(x;a)  -  G^(a)u^(x)|  <  a^lu^(x; 


(3.5) 


< 


a^(a) 


a  (a) 


M 


.  I  o  (a  )u  (a) 


<  e  la  (a)u  (a)  I  , 


I.e.  , 


c^o(a)Uo(x)l  <  jTi"  |u(x;a)|  <  (1  +  2e^)|u(x;a)|  , 


from  which  the  desired  result  follows 


u(x;a)  -  a^(a)u^(x)|  <  e^(l  +  2£^)  |  u  (x;a )  |.  . 


For  k  >  0,  the  same  result  may  be  obtained  by  using  the 
Inequality  (3-5) • 


k 


u(x;a)  -  ^_   a^(a)u^(x) 


< 


n=:0 


'^k+l^^ 


<7v(a) 


k 


a  u 
o  o ' 


k+1 


(x) 


"^(^) 


<  (l+2sjM^^^ 


^k+l(«) 


^i.(a) 


k 


u (x;a) 


<  ej^|u(x;a) 


(3.6; 


In  the  foregoing  discussion  limiting  processes  have  been 
carried  out  in  a  loose  manner  to  save  writing.   The  precise 
meaning  of  inequality  (3-6)  is  given  as  follows:   if  a  (a) 


n 


and  u^(x),  n  =  0,  1,  ...,  N,  satisfy  equation  (3.3)  and  ij)  A) , 
then  for  a  given  e,  >  0,  there  exists  a  a,  >  0,  such  that  for 
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a  t    (0,a,  ) , 


k 

|u(x;a)  -  ^  c7^(a)u^(x)|  <  E^|u(x;a)|  0.6) 

n=o 

k 
for  k  =  0,  1,...,  N.   That  Is,  ^  o    (a)u  (x)  is  an  approximate 

n=o 
solution  to  equation  (3-1)  within  a  relative  error  e,  . 

In  the  sequel,  discussions  will  be  restricted  to  uniformly 
asymptotic  expansions  and  approximate  solutions  within  a  rela- 
tive error. 

However,  uniformly  asymptotic  expansions  are  most  useful 
In  the  cases  where  the  solution  u(x;a)  to  equation  (5-1)  Is  not 
known  a  priori.   In  such  cases.  It  Is  necessary  to  obtain 
Information  about  o    (a)  and  u  (x)  directly  from  equation  (3-1) • 
A  theory  which  can  achieve  this  task  Is  called  a  uniformly 
valid  asymptotic  theory.   In  general.  If  a  uniformly  valid 
asymptotic  theory  Is  realizable,  the  asymptotic  ordering 
a  (a)  will  manifest  Itself  after  certain  properly  arranged 
limiting  processes  are  taken  on  equation  (3-1)  and  the  function 
u  (x)  will  be  governed  by  a  sequence  of  equations  resulted 
from  these  limiting  processes  which  have  the  form 

L  (x,  u  (x) )  =  g  (x)  , 

which  may  be  called  the  n-th  order  asymptotic  equation  for 
the  problem. 
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III. 2.   Nearly  Free -molecule  Flows 

As  stated  in  the  first  chapter,  the  present  study  is 
primarily  concerned  with  the  problem  of  a  uniform  rarefied 
gas  flow  past  a  body  under  the  condition  that  the  character- 
istic dimension  of  the  body  is  much  smaller  than  the  mean  free 
path  of  the  gas.   This  condition  is  called  the  nearly  free- 
molecule  condition.   Throughout  the  analysis  it  is  assumed 
that  the  gas  is  of  infinite  extent  and  the  body  is  a  two- 
dimensional  cylinder  of  convex  cross  section.   Three-dimen- 
sional problems  can  be  treated  in  a  similar  manner. 

Let  the  number  density,  the  temperature,  and  the  velocity 

of  the  undisturbed  flow  be  n  ,  T  ,  and  U/2RT~,  and  the 

o   o      — '    o 

coordinate  system  be  so  chosen  that  the  cross-sectional 
contour  of  the  cylinder  is  defined  by  a  single-valued  and 
piecewise  dif ferentiable  function  R  {^)    through  the  relation 


=  (HiQ) 


This  is  the  very  problem  considered  in  section  II. 5- 
The  mathematical  problem  is  then  to  find  a  function  $(c_;r,05a) 
satisfying  equations  (2.l6)  and  (2.17).   For  this  problem, 
the  nearly  free-molecule  condition  is  simply  a  «  1.   The  task 
is  to  find  uniformly  asymptotic  expansions  for  ^(c_;r,9;a),  i.e., 
to  derive  a  uniformly  valid  asymptotic  theory  under  the  condi- 
tion a  «  1.   By  uniformly  asymptotic  is  meant  uniformly  in 
(r,0)  plane. 
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I 

Before  going  Into  the  details  of  how  such  a  theory  Is 
worked  out.  It  Is  found  clarifying  to  point  out  two  basic 
ideas  which  underlie  the  entire  theory.   For  convenience  of       > 
discussion,  let  the  (r,0)  plane  be  divided  into  two  regions,      , 

an  inner  region  for  r  <  t  and  an  outer  region  for  r  >  t.   As      i 

I 
a  function  of  a ,  t  has  the  property  that  .  . 


11m  T  =  00  (5.7a) 

a  ->  o 


and 


lim  aT  =  0  (3.7b) 

a  ^  o 

while  the  explicit  form  of  T(a)  is  yet  to  be  determined.   The 

two  basic  ideas  are  as  follows.   Firstly,  at  a  point  in  the 

inner  region,  the  distance  from  the  body  is  very  small  compared 

with  the  mean  free  path  l/a ,  as  is  evidenced  by  (3 •7b),  and 

consequently  the  collisional  effects  are  only  secondary,  at 

least  for  molecules  of  moderate  speed.   Secondly,  in  the  outer 

region  the  angle  subtended  by  the  body  is  very  small,  and  in 

fact  if  the  mean  free  path  l/a  is  used  as  the  length  scale, 

in  the  outer  flow  field  the  only  quantities  which  involve  a 

are  the  sizes  of  the  body  and  the  inner  region  which  both 

shrink  to  zero  in  the  limit  a  — >  0.   These  different  features 

in  the  inner  and  outer  regions  seem  to  suggest  that  ^{c_;r  ,6;cl)  , 

in  the  limit  a  — >  0,  has  different  behavior  in  the  two  regions. 

Therefore,  to  investigate  the  nature  of  ^  in  the  limit  a  — >  0, 

it  is  expedient  to  consider  ^(c_;r,0;a)  as  defined  by  two 
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functions:   ^(c_;r,0;a)  restricted  to  the  inner  region  and 
^{c ;p ,6;a)    restricted  to  the  outer  region  where  the  Independent 
variable  r  has  been  changed  to  p  =  ar  as  is  suggested  by  the  use 
of  the  mean  free  path  l/a  as  the  length  scale.   The  precise 
definitions  of  these  new  functions  are 

A 

^(c_;r,0;a)  =  $(c_;r,0;a)   for  r  <  t  (3-8a) 

5(c_;p,0;a)  =  i(c_;  ^,    9;   a)      for  p  >  ar .       (5.8b) 

So  far  no  elements  of  arbitrariness  have  been  introduced 
into  the  above  definitions;  they  are  simply  different  ways  of 
rewriting  ^(G_;r,0;a).   But,  as  mentioned  earlier,  if,  in  the 
limit  a  — >  0,  ^(c_;r,0;a)  has  different  behavior  in  the  outer 
and  inner  regions,  ^(c_;r,0;a)  and  ^(c_;p,0;a)  will  eventually 
have  to  be  treated  as  different  functions.   The  aim  of  the 
uniformly  valid  asymptotic  theory  is  to  obtain  uniformly 
asymptotic  expansions  for  ^(c_;r,0;a)  and  ^(c_;p,0,a),  in  the 
limit  a  — >  0,  directly  from  equations  (2.l6)  and  (2.17).   To 
this  end,  equation  (2.l6)  is  rewritten  in  terms  of  the  new 
functions  ^(c_;r,0;a)  and  ^(c_;p,0;a).   In  the  inner  region 
equation  (2.l6)  then  becomes 

-c         ^ 
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+  0(1,  [fi^^^'J  ^^r^[(^^^^:]  (3.9) 


'^b  rSw 


X  oo 


+  ^rj  J  ir.dr'.^jiT.jf'] 


rsr^vi^ji        ^T, 


where  p'  =  ar '  and  in  the  Integrands  the  following  short -hand 
notations  have  been  used: 


T, 


12 


=  \z(^'}   r,rj,)    =  T,.i(c,;ri';^j 


And,  in  the  outer  region  p  >  aT,  equation  (2.l6)  may  be 
written  as 
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e>o 


where 


-f-^ 


J.[|^T<iK'4-5[-[i'i;df'] 


Rb  «-^ 


C^C 


4-'-    -^       -oo       ^        L«,  -t,  a^     (S)^       -I 


OO  OO       ^    oo       ~<S',  tfxi        <3 


X 


0,-^^m<S),      ,      a^=/5rn@,  , 


and  in  the  Integrands 


A       A 


And,  condition  (2.7)  simply  goes  over  to 

(5.11) 
Thus  far  no  approximation  has  been  Introduced  Into  equation 


11m  f(c_;p,0;a)  =  0    for    c   <  0 
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{^■9),    (3-10)  and  (5-11);  they  are  exact.   But  from  now  on  the 
validity  of  equation  (3-9)  will  be  restricted  to  the  Inner 
region  r  <  t,  and  that  of  equations  (3.10)  and  (3.11)  to  the 
outer  region  p  >  aT  so  that  different  simplifications  can  be 
made  accordingly  to  derive  asymptotic  expansions  for  $  and  J. 


III. 3-   Outer  Solutions 

A 

Let  it  be  assumed  that  for  the  functions  ^(c_;r,0;a)  and 
^(c;p,0;a)  there  exist  uniformly  asymptotic  expansions  of  the 
type  Zo^{a)'^^{c_;r,9)    and  E  a"^(a)_^^(c_;p ,  e)  which  satisfy 
conditions  (3-3)  and  {^A).      It  will  be  shown  that  under  this 
assumption  the  fun  tlons  a  (a),  a  (a)  and  J  (c;r,0)  can  be 
derived  from    the     equations  (3-9)  and  (3.10).   Also 
derived  from  (3-10)  is  a  sequence  of  integral  equations,  each 
Independent  of  O^,  whose  solutions  are  the  functions^  (c_:p,e). 

Since  both  ^  and  ^  appear  in  equations  (3-9)  and  (3-10) 
it  is  necessary  to  make  certain  estimates  for  one  or  the  other 
to  initiate  the  asymptotic  analysis.   The  estimate  which  is  made 

A 

here  is  that  the  lowest  order  asymptotic  expansion  of  ^  in  the 
inner  region  is  the  free-molecule  solution  to  the  problem. 
With  this  estimate  0^^(01)  can  be  explicitly  obtained  from 
equation  (3. 10),  which,  in  the  limit  a  — >  0,  becomes  an  equation 
for  ^^  and  is  free  of  a.   Making  use  of  these  results  for 
a^(a)  and  ^^,  ^0^0  *^^"  then  be  evaluated  from  equation  (3-9) 
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and  Is  to  be  checked  against  the  estimate  previously  made  to 
ensure  self -consistency .  In  what  follows,  the  results  for  o 
^  will  be  worked  out  In  steps. 

First  of  all  the  asymptotic  expansions  of  the  known  func- 
tlon  ^   (c_;  ^,    9 ,a)    and  the  various  integrals  involving  a  $ 
are  to  be  calculated  in  the  outer  region.   In  this  connection, 
it  is  noted  that  in  the  outer  region  Qj)-,  and  ®^   are  very  small 
angles  and  can  be  expanded  asymptotically  as 

r 

®2(&'  «)-F  ®20('"  -^^  ®21«"    ■ 

It  will  then  be  shown  that  the  function  J*  (c_;  ^,  0;   a) 
can  be  expanded  asymptotically  as 

^  (c;  ^,  0;  a)^  f  %{^;P  ,6)     • 
To  this  end,  it  is  observed  that 

q;i^,i  ,»j'^)  -■.  -2uI,[cos&„{c,+fU^:i)]  (3.15) 

In  the  above  integrand  the  factor  cos  6      calls  for 
special  precautions  because  of  its  strong  dependence  on  't>    . 
To  overcome  the  difficulty  caused  by  this  strong  dependence, 
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it  Is  useful  to  Introduce  a  new  coordinate  system  {^,0)    for 

the  cross-sectional  contour,  where  0  Is  the  polar  angle  0   of     | 

{p ,6)    and  p  Is  the  angle  between  the  position  vector  of  a        ( 

point  (p,0')  on  the  contour  and  a  straight  line  of  Inclination 

6.       (See  Figure  5-)   In  the  {^,0)    coordinates,  R,  ,  0'  and  6 

may  be  written  as  p  =  p(p,0),  0'  =   0'((3,0),  and  0=0  (p,0). 

i 
In  the  limit  of  t  — >  oo  while  p/a  >  t,  a  simple  form  can  be       ' 

found  for  *   as  a  function  of  p  with  fixed  0  and  p.   In  fact,     | 

.    ap (p, 0)   .   „ 
1      p        I-    } 

! 

^\   =  p  c'^s  0^  (P  ^°^  P  +  sm  p  §)dp  ,  , 


a    1 


p  cos  4)-, 
and  as  p/ct  — >  <», 


6(P,0)dp  , 


*1^  p  P(P'^)sln  p  , 

d*^^^  (p  cos  p  +  sin  p  ^)dp 


^O,0)dp 


For  given  (p,0)  the  values  of  p  which  correspond  to  <t>-.    equal 
to  -  @-,  and  (h)     are  respectively 

P^(p,0)=  sln-^^^i-^)^-  sm  @)^)  (3.14a) 


and 
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^2^P>e)   =   sln-l(^^^sln@2)  (3.l4b) 

Then,  upon  changing  the  Integration  variable  0-,  Into  p, 
equation  (3-13)  becomes 

i^'J.^j^)  ^  -  f^  J  ^/3  r[^.fl^r;(/i^^^ -Jfc^sA)}!       (3.15) 

In  the  limit  p^'a  — >  <»,  P-,(r,0)  and  Pp(r,0)  may  be  represented  by 
p^(r,e)::^sln-l(^^(H)^^(0)) 


or,  on  solving  the  above  equations. 


p^(r,0)^  sln-^(^  ©10^^^)  (3.16a) 

P2(r,0)^  sln-^(^  (B^q{9))  (3.l6b) 


which  will  be  written  respectively  as  ^-.(S)  and  Pp(0). 

If,  in  obtaining  the  asymptotic  expansions  for  ^  {c;   ^,  G;   a) 
the  limits  of  integration  -P-,(r,0)  and  f52i^>^)    i"  equation  (3-15) 
can  be  replaced  by  -^-.{0)    and  pp(e)  as  defined  by  (3.16),  then 
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^(^-,h;^j  =  -'^ Af  ru^t^-i i^f¥^/  --f'-'^J}]    \ 


^  'P,w 


In  the  limit  a  ->  0,  the  above  Integral  is  found  to  approach 
a  function  Independent  of  a .   If  this  function  is  denoted  by 
p- Ip  {c_;p  ,9;a) ,    the  lowest  order  asymptotic  expansion  for 
^  (c_;  ^,0;a)   may  be  written  as 

^  (c;  £-,  0;  a)^-  ^^^(c;p,e)  .  (3.1?) 

The  derivation  of  (3-17)  is  lengthy  and  may  be  found  in 
Appendix  II.   Asymptotic  expansions  of  this  type  are  concep- 
tually analogous  to  the  multipole  expansions  for  a  localized 
electrostatic  charge  distribution,  the  lowest  order  expansion 
being  the  monopole  term  or  the  contribution  due  to  the  net 
charge.   Following  the  same  method  used  in  Appendix  II,  the 
replacement  of  the  limits  of  integration  -P-,(p,0)  and  p  (p,0) 
in  equation  (5-15)  by  -^-,(0)  and  ^2^^^'  defined  in  (5-16), 
can  be  shown  to  be  legitimate. 


In  a  similar  manner,  the  lowest  order  asymptotic  expan- 

sions 

1.-  p.  ~  "      ■  - 

obtained  as 


for  B^[  i  ^Tt,Jr'^^  f  ^Ti^^p'l  ,     i  =  1,  2,  and  I^  are 


/?b         ^c     J  J  (3.18) 

and 
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X|  OJ  ^  J^  1,0  Ij  1 


where 


R^  =  R^O,p,0) 


00  >0 


^,(^) 


It  should  be  noted  that  the  asymptotic  expansions  (3-17) 
and  (3.18)  are  only  valid  for  c  not  too  large.   More  precisely, 
they  are  only  valid  for  c  <  Ap(a)  where  the  function  A2(a) 
approaches  Infinity  as  a  — >  0  (cf.  Appendix  II).   This  restric- 
tion is  not  detrimental  to  the  asymptotic  theory,  for  there  is 
only  an  exponentially  small  number  of  molecules  with  a  velocity 
c  >  Ap(a)  and  consequently  the  contribution  from  such  molecules 
to  any  moments  of  the  distribution  function  is  exponentially 
small  and  thus  altogether  negligible.   To  obtain  higher  order 
asymptotic  expansions  f or  ^  (c_;r,0;a),  B-^,    and  B^,    it  is 
necessary  to  take  into  account  the  effects  due  to  the  persis- 
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tence  of  the  molecular  velocity  and  the  much  shorter  free  path 
of  the  slow  molecules;  neither  of  these  effects  are  appreciable 
In  the  lowest  order  expansions  (3-17)  and  (3-18). 

Now  the  lowest  order  inhomogeneous  term  In  equation  (3-10) 
Is  of  the  order  a  and  thus  suggests  that  cj  (a)  -   a.   Let 
equation  (3.10)  be  divided  by  a .   If  a  limiting  process  a  — >  0 
is  performed  on  this  new  equation  and  if  the  order  of  this 
limiting  process  and  the  integrations  can  be  interchanged 
freely,  the  following  equation  will  result: 

^  ^    r      ^      Ki,  (3.18) 


where 


In  the  above  expressions,  for  all  the  terms  involving 
a,  T,  ®   and  @p,  a  limiting  process  a  — >  0  is  understood 


to  have  been  taken.   To  carry  out  this  process  more  explicitly, 
it  is  necessary  to  know  the  behavior  of  5  as  p  — >  0.   How- 
ever,  since  a  $  and  a  ^  must  match  asymptotically,  it  is 
seen  that 

I  ^   ^.(£>P^Q)   as   p  ->0  . 
^o       p  ^ 
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Now  Hadamard ' s  finite  part  concept  can  be  used  to  carry 

out  the  limiting  process.   For  example,  the  operands  of  the 

operators  B   ,  in  the  limit  a  — >  0,  may  be  calculated  in  the 
following  manner 

/'^(i,T,cir'  ^l^S%M'j''^^^i-^)\^^'       (5.19) 


A 


where 


and   oo 

A   A 


which  is  a  finite  quantity,  independent  of  a. 
And  similarly 


oo 


where 


and  <^^  ^ 


(:5.20) 
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Combining  equations  (3-19)  and  (3-20)  yields  the  following 


results 


^b 


'A  ^-C  ' 


Therefore,  as  a  — >  0,  the  terms  of  B.q  and  B^q  approaches 
zero;  i.e.. 


Oo 


In  a  similar  manner,  it  may  be  shown  that 
a— ^  o  ^^ /?b 


oiX^ 


and 


lim  Ii^[$]  =  0 


a  ^►o 

Now,  consider  the  term  in  (3-l8)  which  Involved  I^.. 
the  variable  *   is  changed  into  p  through  the  relation 

p  sin  p  =  p ' sin  0   , 


If 


this  term  may  be  written  as 


P 


0    "% 


"i  o 


O     6> 
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Finally,  it  can  be  shown  that 


^ztjij.^?']-    l2of/CT,df'] 


where 


-Do    o    L--^   0   -" 
Summing  all  these  results  for  equation  (3.l8)  yields 

the  lowest  order  asyinptotlc  equation  for  ^  as  follows 


/■x^ 


iXD 


,(^.Ae;=-f  if/£.-f.ej  +J3[f  r^TJ/'J  (3.21) 


0= 


-f- 


iJlfoWj^-i^.Lfl] 


/ 


The  boundary  condition  (3-ll)j  when  applied  to  ^  ,  becomes 


lim  f  (g_;p,0)  =  0    for    c  <  0  . 


(3.22) 


Equations  (5-21)  and  (5.22)  define  the  lowest  order  outer 
solution  ^„(c_;p,0).   It  should  be  emphasized  that  equation  (3.21) 
Is  Independent  of  a. 

Also,  It  should  be  noted  that  the  outer  region  Is  not  a 
region  throughout  which  collisional  effects  dominate.   In  fact, 
in  the  part  of  this  region  where  p  — >aT,  which  is  very  small, 
the  nature  of  the  flow  field  is  that  of  the  free -molecule  flow. 
In  the  part  of  the  outer  region  where  p  —^m,    the  flow  field 
is  presumably  a  collision-dominated  one.   And,  in  between,  a 
full  transition  may  be  expected. 
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III. 4.   Inner  Solutions 

After  the  results  have  been  obtained  for  'a      and  ^, 
equation  (3-9)  can  be  used  to  calculate  the  inner  solution 
a^^-,.   First  of  all,  the  inhomogeneous  term  *J^  (G_;r,6,a)  is 
expanded  asymptotically  for  small  ar  as 


.  ^  (c;r,0;a)^  J/^(c;r,e)  +  a\{c_;v ,Q) 


(3.23) 


for  1  <  r  <  T,  where 


^(c;r,e)«.  ^^^[o_;0,e)      as   r 


as   r  -^   00 


5^(c;r,0)<^  G-^(g;0) 

Now,  on  the  R.H.S.  of  equation  (3-9) j  ^  and  J  may  be 
replaced  by  ^  and  a  ¥  respectively.   For  example,  consider 


the  term  I^. 


h^ 


T 

r    ^ 


I  T-j^dr  '  +  -   I  T-]_dp  '  ] 


(5.24) 


ax 


h^ 


r- 


^  T-,  dr  '  +   J^T-j^dp  '  ] 


ax 


which  is  a  consequence  of  the  basic  assumptions  that 

A  A  A 

If    -    $ol     <     ^nl^rl  fO^       r    <    T, 


o'^o' 


$   -  0^  1^1    <  ote    1^    I      for  p   >  ax 
J-  O"*©  '  o '  -'-o '  ^ 
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and  e„(ct)  — >  0  as  a  — >  0. 

For  r  <  T,  the  function  T.  (c_;r  ,r '  ;a) ,  1  =  1,2,  has  the 
asymptotic  expansion 


T.  ^  T 


r 


/r   -r  sm  4). 


(3.25a) 


And  for  T.(c-,;r,  ^  ;a)  It  will  be  assumed  that  ar  «  p', 
and  the  result  Is 


Ti^^l'^'  5-'^) 


T  =  exp{-^p.j 


(3.25b) 


It  should  be  noted  that  (5.25b)  is  not  valid  for  r  -^  i   and 
p'  — >aT.   Therefore,  to  use  (5.25b),  It  Is  necessary  to  make 
the  restriction  that  r  <  t-,  where  t,  has  the  properties  that 


00,  aT-|  — >  0,  and  t-,/t  — >0  as  a  — >  0.   Since  \     d 


■1 


n 


oes 


not  depend  on  a,  the  asymptotic  expansions  2o  ^  for  ^  will 
be  extended  to  cover  the  Interval  (t-,,t). 

With  the  results  in  (3.25a)  and  (3.25b)  for  the  function 
T. ,  equation  (3.24)  may  be  written  as 


r  -  1  ^ 


I  [   ^T^dr'  +^    $T^dp']«iI  [   f  ^T  dr '  +    f^Tdp'] 

•^     ■J  -J  ■^     J  ^J 


aT 


QT 


Now  the  concept  of  Hadamard ' s  finite  part  will  be  used 
to  carry  out  the  limit  a  — >  0  for  the  last  equation.   The 
result  is 


-  ^5 


o<-»o  ^  ^    ^      ^  ^^ 


y^     1,[    iijJr-^fijJf]  (3.26) 


where 


OO 


The  derivation  of  equation  (3-26)  runs  parallel  to  that  for 
equation  (3-19)  In  the  last  section. 

For  the  rest  of  the  terms  on  the  R.H-S.  of  equation 
(3.9),  similar  results  are  obtained  as  follows. 


r 


A  A 


C>o 


00  CO 


(3.27a) 


(3.27b) 


(3.27c) 


-  ^2  C  f  iT^'''liT^f']-(^^JrjQ^cc^ajJ 


00 
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where 


OO       oo   (5), 


o 

2    2 

S-,  C-,    c 

o 


and  R,  Is  a  function  of  (0  ,r,e)  defined  by  the  solution  to 
the  simultaneous  equations 


R. 


^(«b) 


-1  r 
+  *-,  -  sin   (^  sin  <t>  ) 

b 


Substituting  equations  (3-26)  and  (3-27)  into  (3-9) 
yields,  in  the  limit  a  — >  0,  the  following  equations 


A 

0-1 

o 


$   =  $ 


0-,    -    -a  ina    =  ain(— ) 
1  a 


(3.28a) 
(3.28b) 
(3.29a) 


A  .  ^ 


^1  =  V^o^^o'^)]  +  W^o^^'")] 


+  l^[G^{^;e)]   +   l^[G^{o_^,9)]    , 


(3.29b) 


Op  =  a 


(3.30a) 


^7  - 


A  A 


(3.30b) 


-f- 


Equations  (3-28)  are  the  result  of  free-molecule  solution 
and  thus  ensure  the  self -consistency  of  the  theory.   Equations 
(3-29)  and  (3-30)  give  the  first  and  second  order  inner 
solutions . 

After  these  functions  have  been  obtained,  it  is  then 
possible  to  go  to  the  next  order  outer  solution,  which  in 
turn  can  be  used  to  obtain  the  next  order  inner  solution.   By 
continuing  this  procedure,  the  inner  and  outer  solutions  can 
be  obtained,  at  least  in  principle,  to  all  orders.   In  this 
connection,  most  of  the  difficulties  are  due  to  higher  order 
asymptotic  expansions  for  known  functions.   Because  of  diffi- 
culties of  this  sort,  to  find  the  inner  and  outer  solutions 
to  orders  higher  than  what  have  been  obtained,  it  becomes 
necessary  to  divide  the  velocity  space  into  several  regions 
and  to  treat  each  region  properly.   The  difficulties  are 
mainly  caused  by  the  persistence  of  molecular  velocity  and  the 
slow  molecules.   However,  no  attempt  is  made  here  to  carry  out 
the  higher  order  analysis. 
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A 

The  zeroth  order  Inner  solution  ^  agrees  with  the  zeroth 
order  outer  solution  In  a  region  where  r  »  1  and  ar  «  1. 
This  region  may  be  called  the  overlap  region.   From  the  derlva- 
tlon  of  (5-29)  It  Is  clear  that  ^-,  Is  produced  by  the  collisions 
In  the  overlap  region.   The  fact  that  cj-,|-.  Is  of  a  higher  order 
than  o-„£p  may  be  attributed  to  two  features  of  the  overlap 
region,  namely  (l)  the  overlap  region  Is  of  considerable 
extent,  and  (11)  molecules  produced  by  collisions  In  the  over- 
lap region  have  little  chance  of  suffering  a  collision  when 
they  travel  to  a  point  In  the  Inner  region. 


III. 5-   Aerodynamic  Forces 

With  the  results  obtained  In  the  last  two  sections, 
asymptotic  expansions  for  the  distribution  function  f  can  be 
easily  calculated  from  equations  (2.15)  and  macroscopic 
quantities  of  aerodynamic  Interest  can  be  obtained  by  Inte- 
grations over  the  velocity  space.   For  example,  the  normal 
pressure  at  a  point  (a,0)  on  the  surface  of  the  body  (see 
Pig.  •4-)  iTiay  be  computed  according  to  the  following  formula: 

K,  -  -— -vT- JJj  e     cUh^)<^^C  (3. 31) 
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And  It  will  be  assumed  without  loss  of  generality  that  the 
velocity  space  at  {a.,0)    only  consists  of  region  I  and  region 
II,  such  as  In  the  case  of  a  circular  cylinder  In  a  uniform 
flow.   For  more  general  body  shapes,  the  presence  of  region 
III  causes  no  essential  difficulties.   For  a  circular  cylinder, 
a  =  1,  and  equation  (5-31)  may  be  written  as 


P.fiTo  TT 


c-„>^ 


^    ^ ' '     ('€') c^[^~ ^. ^.  4  ( ^. 'S-^JJ  ^^c 


^  fi^JJl^'P 


C„<0 


If  the  inner  and  outer  solutions  are  substituted  into 
(3-15a.)  and  (3-15b),  the  last  equation  becomes 


Po/fr     '         ^■^'■^ 


/e"  'e 


A      A 


Again  by  using  the  concept  of  finite  part,  the  limiting 
process  a  — >  0  can  be  carried  out  for  equation  (3-32).   The 
result  is 


3. 


-\^-^(^^0^^'^-  (3.33) 
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similar  results  can  be  obtained  for  the  tangential  stress 
P  at  (1,0) .   In  fact, 

^^  -/  =  -  ^^;n6^  .  ^JLlk)i^^^  I  ^„fte;>     (3.54) 


A   A 


+  ^<c„^,     (LTd.Vf^rj/'  > 


with  these  results,  the  drag  can  be  computed  according 
to  the  formula 
27r 


D  =  - 


(P  cos  0  +  P, sin  9  )d( 
^  n     n    t     n'^ 


o 


For  a  circular  cylinder 


n 


',  and  the  expression  for 


D  can  be  written  more  explicitly  as 


D   =  ^^^  +  a    -^n(|)Di   +  cxD^    , 


where 


°fM  ^    ^2^   (||+    1)   +    1]U 


^^^  -^\      (w(c;e)|    GQ(c;0))de 


StT  00  00 
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and 

wfc:0l  =       ^ .  ,,    „  . 

./—  n  t        ^4    4   ,/—  n 
yir  yw 


2 
w(c_;0)  =  ^3  c^c^sln  0  +  {-^ if  +  "z:  c^)cos  0  . 


Again  It  Is  noted  that  the  first  order  correction  term  D-, 
is  due  to  the  collisional  effects  in  the  overlap  region  where 
the  inner  and  outer  solutions  match  asymptotically.   It  is  also 
noted  that  D-,  can  be  calculated  by  Willis'  integral  iteration 
method.   For  practical  purposes,  the  number  ^n(l/a)  may  not  be 
very  large  and  therefore  for  reasonable  accuracy  it  may  need 
to  take  the  combination  a  ^n(l/a)D-,  +aDp,  instead  of 
ain(l/a)D-|,  to  account  for  the  collisional  effects.   In  such 
circumstances,  ^  has  to  be  computed  numerically  to  give  the 
result.   All  that  has  been  said  in  this  paragraph  regarding  the 
drag  D  is,  of  course,  true  for  general  two-dimensional  flows  and 
is  not  restricted  to  flows  past  a  circular  cylinder. 


III. 6.   Question  of  the  Mathematical  Foundation 

It  is  perhaps  clarifying  to  point  out  here  the  logical 
implications  of  what  has  been  done  in  the  last  two  sections. 
The  results  therein  may  be  summarized  in  the  following  statement. 

Statement  A.   If,  for  (^  and  ^  there  exist  asymptotic  expan- 
sions  of  the  type  Za^(a)  $  (c_;r,e)  and  Zo    (a)  ^  (c_;p,0) 
satisfying  conditions  (3-3)  and  (3-^)  with  the  property  that 
£y$Q  and  $V^„  are  uniformly  bounded  in  r  <  t  and  p  >  ax 
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respectively,  and  If  o  ^  Is  the  free-molecule  solution,  then 
the  functions  o    ,   'o    ,    '^    ,    and  ^   can  be  obtained  by  the  procedures 
described  in  sections  III  3  and  III. 4,  and  in  particular,  for 
two-dimensional  problems, 

(i)  a   =  1,  a   =a^n(l/a),  and  a^   =  a 
(ii)  ^-1  =  R.H.S.  of  equation  (3-29b) 
(iii)  ^2  =  f^-H-S.  of  equation  (3.30b),  and 
(iv)  f   is  the  solution  to  equations  (3-21)  and  (3.22). 

There  are  two  conditions   in  Statement  A,  namely  (a)  for 

^  and  ^   there  exist  uniformly  asymptotic  expansions  of  the  type 

Zo   $  and  2  a  (b   satisfying  conditions  (3-3)  and  (3-'^),and 
n^n        n^n       jo 

A 

(b)   a  $   is  the  free-molecule  solution  to  the  problem.   Prom 

a  rigorous  viewpoint,  these  two  conditions,  if  true,  must  be 

proved  as  a  consequence  of  the  original  governing  equations 

(2.l6)  and  (2.17).   In  most  asymptotic  theories,  this  is  a  very 

difficult  task  and  in  the  absence  of  a  proof  the  conditions  (a) 

and  (b),  in  one  form  or  another,  are  usually  taken  together 

with  the  existence  and  uniqueness  of  the  solution  as  the  basic 

assumptions  of  the  theories.   One  way  to  prove  conditions  (a) 

and  (b)  is  to  verify  the  converse  of  Statem.ent  A  which  may  be 

called  the  Asymptoticity  Theorem  and  is  stated  as  follows. 

Asymptoticity  Theorem.   If  the  functions  a  ,  a"  ,  $  ,  and 

^  are  obtained  by  the  procedures  described  in  Sections  III. 3 

/\  — < 

and  III. 4,  then  Sa  ^  and  Sa  ^  are  respectively  the  uniformly 

asymptotic  expansions  for  $  and  ^  satisfying  conditions  (3-3) 

and  (3-4) • 
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This  theorem.  If  proved  true,  together  with  the  existence 
and  uniqueness  of  the  solution  to  the  problem  will  be  sufficient 
to  furnish  a  mathematical  foundation  for  the  asymptotic  theory. 

Since  a  proof  of  the  Asymptotlclty  Theorey  has  not  been        ; 

i 

bound  for  the  general  boundary  value  problems  for  the  (linearized) 
Boltzmann  equation.  It  appears  desirable  to  apply  the  proposed 
asjnnptotic  theory  to  a  simple  problem  for  which  the  Asymptotlclty 
Theorem,  as  well  as  the  uniqueness  and  existence  of  the  solution, 
may  be  rigorously  proved.   It  Is  hoped  that  the  application  to 
this  simple  problem  will  help  to  illustrate  the  mathematical 
nature  of  the  theory.   The  simple  problem  which  is  to  be  examined 
in  the  next  two  chapters  is  that  of  a  circular  cylinder  rotating 
in  an  infinite  gas  which  is  at  rest  at  infinity.   And,  since  the 
proposed  asymptotic  theory  is  not  restricted  to  the  Boltzmann 
equation,  a  simpler  model  equation  will  be  used  for  the  cylinder 
problem  for  the  sake  of  further  simplicity. 
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rv.   The  Cylinder  Problem--Asymptotic  Theory 
rv.l.   The  BGKW  Model  Equation 

The  problem  to  be  studied  In  this  and  the  next  chapters 
is  that  of  a  circular  cylinder  rotating  in  an  infinite  gas 
which  is  at  rest  at  infinity.   For  this  problem  a  model  equa- 
tion proposed  independently  by  Bhatnagar  et  al .  (195^)  and 
Welander  (195''^)  is  used  instead  of  the  Boltzmann  equation. 
This  model  equation,  which  may  be  called  the  BGKW  equation, 
embodies  the  assumption  that  the  rate  of  change  of  the  distri- 
bution function  F  due  to  collisions  is  proportional  to  F  -F, 
the  deviation  of  F  from  the  local  equilibrium  value;  that  is 


/coll 


where 


1   5/2     f   (i-q)^l 
^e  =  "(2¥rt)     ^^P\-  -2K^j 

and  q,  T,  and  n  are  functions  of  position  only,  defined  by 
the  following  equations 

n  =  ;//  Fd^^  , 

and       nRT  =  ^  ///(i-q)^  Fd^^ 

and  the  proportionality  coefficient  a  is  in  general  also  a 
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function  of  position.   The  BGKW  equation  can  then  be  written  as 

^=a{F^  -P)  .  (4.1) 

It  is  seen  from  the  definitions  of  q,  T,  and  n  that  the 
conservation  laws  of  matter,  momentum  and  energy  are  main- 
tained collectively  by  the  collisions.   Taking  a  as  the  local 
collisional  frequency,  the  BKGW  equation  may  be  derived  from 
the  physical  argument  that  the  molecules  emerging  from  a  col- 
llcion  adjust  themselves  instantaneously  to  the  local  equilib- 
rium state  of  the  gas  because  collisions  in  general  randomize 
the  distribution  of  velocities  among  the  colliding  molecules. 
Although  this  argument  sounds  physically  plausible  a  justi- 
fication of  the  use  of  the  BGKW  equation  as  an  approximation 
to  the  Boltzmann  equation  for  some  rational  molecular  model 
is  difficult  to  arrive  at  from  a  rigorous  viewpoint.   However, 
since  this  equation  is  used  here  only  to  help  illustrating 
the  nature  of  the  asymptotic  theory  proposed  in  the  last 
chapter,  a  different  viewpoint  may  be  taken  toward  the  question 
of  its  justification.   In  this  connection,  the  following  known 
results  based  on  this  equation  may  be  cited.   The  Ghapman- 
Enskog  expansion  scheme  has  been  applied  to  this  equation  for 
a  »  1,  e.g.,  in  Llepmann  et  al .  (I962),  and  to  the  zeroth 
order  of  a~  yields  the  Euler  equations,  which  is  Identical  to 
the  results  for  the  Boltzmann  equation.   The  first  order 
results  are  very  similar  to  those  for  the  Boltzmann  equation 
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although  they  give  an  Incorrect  value  of  the  Prandtl  number. 
For  some  low-speed  linearized  gas  flows  the  temperature  field 
is  not  coupled  with  the  velocity  field  and  the  incorrect 
value  of  the  Prandtl  number  may  not  be  a  serious  defect.   On 
the  other  hand,  for  a  «  1,  the  BGKW  equation  has  been  used 
for  the  linearized  Couette  flow  by  Willis  (1958)  and  shows  that 
the  macroscopic  velocity  gradient  half-way  between  the  plates 
is  proportional  to  ain(l/a).   This  result  is  in  agreement 
with  that  obtained  for  the  Boltzmann  equation  for  hard  spheres, 
cf.  Willis  (1961).   These  are  among  the  few  results  for  which 
a  comparison  between  the  BGKW  equation  and  the  Boltzmann 
equation  for  some  molecular  models  is  available,  and  they 
seem  to  indicate  that  the  BGKW  equation  is  a  qualitatively 
reasonable  approximation  to  the  Boltzmann  equation,  at  least 
for  the  linearized  flow.   This  indication  also  seems  to  have 
been  reinforced  by  the  results  of  Gross  and  Jackson  (1959)- 
The  linearized  Boltzmann  equation  may  be  written  as 


df     r     1 


2 


e   ^  k(c,c^) f  d^c^  -  f ]  ,   (^.2) 


where 

k(c,c^)  =1+1  (c^-  |)(c^  -  |)  +  2  c,c^  .     (4.5) 

The  notation  is  identical  with  that  introduced  in 
Chapter  II.   Equation  (4.2)  may  also  be  written  as  the  fol- 
lowing system  of  equations. 
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df 

dt 


and 


-   a[(7r)-^/2  [-^  +  (c2-  |)^  +  2c. u]  -  f  j    {^  A) 
I  o  o 

2 

n  =  n^  ///  e   ^  f d^c^  , 

2 
T  =|t^  ///  e"""!  (c^  -  |)  fd^c  , 

2 


IV. 2 .  Formulation  of  the  Problem 

For  the  general  boundary  value  problem  for  a  uniform 
flow  past  an  immersed  body,  the  governing  equation  can  be 
obtained  by  substituting  equation  (il-.^)  into  equation  (2.l6). 
However,  for  the  cylinder  problem  It  is  much  simpler  to 
obtain  the  governing  equation  directly  without  using  equa- 
tion (2.l6).   For  this  linearized  problem,  n  and  T  are 
second  order  quantities  in  the  angular  speed  of  rotation, 
as  can  be  seen  by  reversing  the  direction  of  rotation,  and 
may  therefore  be  neglected.   The  BGKW  equation  then  becomes 


^  =   a[2(7r)-^/2  ^.^  _  ^j  (4,5) 


where 


2 

R=  III  ^l   e      ^   f d^c^ 


If  the  origin  of  the  (r,0)  coordinate  system  is  placed 
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at  the  center  of  the  cross  section  of  the  cylinder,  the 
r-component  of  the  mean  velocity  of  the  gas  vanishes  iden- 
tically as  a  consequence  of  the  conservation  of  molecules 
and  the  boundary  conditions  at  the  cylinder  surface  and  at 
infinity.   Therefore  equation  (^-5)  can  further  be  simplified 

to 

2 


r  or    r  oc      r   dc^    l  v  /      y   \  / 


where  a  is  a  constant  inversely  proportional  to  the  mean  free 
path  of  the  gas  and  a'>(r)  is  the  mean  angular  speed  of  the  gas 
defined  by 

C  C  C  2 

a3(r)  =  ^     e"^  f  c^d^c  .  (4.?) 

J  J  -^ 

The  boundary  conditions  for  equation  (4. 7)  are 


and 


f  =  0   as    r  — >oo,   for  c^  <  0  ,  (^-8) 


f  =  2{ir)   ^/^  Cg   at  r  =  1,   for  c^  >  0  .        (4.9) 


In  obtaining  the  boundary  condition  (4.9),  both  the  radius 
and  the  angular  speed  of  the  cylinder  have  been  set  to  1 
for  the  sake  of  simplicity. 

Equation  (4.6)  and  equation  (2.12a)  have  the  same 
differential  operator  for  which  the  characteristic  equations 
are  given  by  (2.13)-   Along  these  characteristic  curves  f  may 
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be  integrated  In  terms  of  a)(r),  using  equation  (^-6)  with 
the  boundary  conditions  (4.8)  and  (4.9).   The  results  are 

For  region  I  (tt/S  <  |*|  <   w) , 


<^l[0^^-'H(f^^ -'^'^m^^  .  (4.iOa) 


for  region  II  (|0|  <  sin~  — )  , 


f  ■=  ^Srs!^l4{^xf[-^rco^<f-  /-r^^/n^^  )]  (4.iob) 


-h 


TT 

r 


TT 


i/Z 


r6ir)ic^l  j'-^=3==r  e,j^j-^(rCo^^-L^,r^S,'^y)  \ 


^  r'6in  y 


-1      1  TT 

and   for   region   III    (sin       —  <    |  *  |    <  y^) , 


r  ^  II  ^2 
oo 


u 


where   s  =  /c  +  Cg  and  4)  =  tan"  cVc  . 

Substituting  equations  (4-.  10)  into  equation  (4.7)  gives 
a  linear  integral  equation  for  a3(r).   From  physical  considera- 
tions the  unknown  function  co(r)  will  be  restricted  to  be 
uniformly  bounded  in  r.   For  such  functions,  the  order  of 
integrations  over  s,  4),  and  t)  in  the  linear  integral  equa- 
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tlons  for  co(r)  may  be  Interchanged  according  to  the  Fublni 
theorem  in  the  Lebesgue  Integration  theory  and  the  result 
may  be  written  as 


CO 

cu(r,a)  =  fi(r,a)  +    K(r,Ti  ;a)aj(ri  ,a)dT] 


(4.ii: 


where  the  dependence  of  a)(r)  on  a  has  been  expressed  explicitly 
and  the  functions  n(r,a.)  and  K(r,ri;a)  are  given  by 


t  ..... 


JJx]  = 


r        2 

.-  -s 
s  e 


n  _  -  s  "  -  x/  s 


ds 


o 


-1  n 
b  =  sin   —    for   ■(]  <   r    , 


b  =  2 


for   11  >  r 


The  boundary  condition  at  Infinity  requires  that 


11m  cD(r,a)  =  0 
r  ->  CO 


(^.12) 
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The  asymptotic  theory  proposed  in  Chapter  III  will  be 
applied  to  equation  (^.11)  under  the  boundary  condition 
(^.12).   In  the  next  two  sections,  the  Inner  and  outer  solu- 
tions will  be  obtained.   That  the  Inner  and  outer  solutions 
are  Indeed  uniformly  valid  asymptotic  solutions  to  the  problem, 
as  well  as  the  existence  and  uniqueness  of  the  exact  solution, 
will  be  proved  In  the  next  chapter. 


rv.3-   Outer  Solutions 

As  proposed  In  Chapter  III,  the  physical  space  Is  divided 
Into  an  Inner  region  for  r  <  t  and  an  outer  region  for  r  >  t, 
where  t  — >  <»  and  qt  — >  0  as  a  — >  0.   In  the  Inner  region 
cD(r,a)  Is  denoted  by  (Xi(r,a),  while  In  the  outer  region  p  =  ar 
and  ^  =  ar;  are  used  respectively  as  the  Independent  variable 
and  the  variable  of  Integration  and  co(^,  a)  Is  denoted  as 
cD(p,a).   The  task  of  the  asymptotic  theory  Is  to  seek  the 
asymptotic  solutions  2  a  (a)a)  (r)   and  Z  'o    (a)co'  (p )  for  a)(r,a) 

and  co(p,aj  under  the  condition  a  «  1. 

„. — 

For  the  function  cD(p,a)  the  governing  equation  (4.11) 

and  condition  (4.12)  become  respectively 


(4.13) 


-  62  - 


oCZ^ 


and 

lim  co(p,a)  =  0  (^-1^) 

p  ^  00 

Now  the  asymptotic  expansions  So  (a)cu  (p)  will  be 
derived  from  (^.13)  and  (4.l4)  for  co'(p,a)  In  the  outer 
region  p  >  ax. 

The  Inhomogeneous  term  fi(— ,  a)    has  the  following  asymp- 
totic expansion,  for  ^  »  1, 

For  the  derivation  of  (4.15),  see  Appendix  III.' 

The  kernels  K(^,  — ;  a)  and  K(^,  r\;a)    can  be  simplified 
by  making  use  of  the  fact  that  ^  >  t  and  t  — >  oo  as  a  — >  0. 
For  example.  In  the  limit  a  — >  O  K(^,  — ,  a)  becomes  aK(p,^) 
where  K(p,^)  Is  Independent  of  a.   The  explicit  form  of 
K(p,e)  Is 


i 


J/ai^^ 


•"  lq^JiJ:^^'^f"'^^W^]^f  \ 


where 


-1  6 
b  =  sin   ^    for    ^  <  p  , 

r 


b  =  |-         for    ^  >  p 
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Equation  (^.l6)  can  readily  be  derived  from  the  expres- 
sion for  K(— ,  — ;  a)    by  carrying  out  the  limiting  process 
a  ->  0. 

For  the  kernel  K(— ,  n ;  a)  It  will  be  assumed  that 
ri  «  £-.   And  consequently.  In  the  limit  a  — >  0, 


where  t-.  =  sin"  ^  and  tp  =  sln~  — ,  or,  more  explicitly, 

/^^^'?  -;  -  y^frr-  IJJ^.  .5.-^  Jirf]       '""' 

It  Is  noted  that  the  condition  n  «  £-  and  therewith 

a 

equation  (^.17)  are  not  valid  for  the  whole  outer  region, 

since  p  may  approach  qt.   Therefore,  to  use  expression  (4.17), 

It  Is  necessary  to  make  the  restriction  that  p  >  aT.,  where 

T-,  (a)  has  the  property  that  t,  — >  oo,  aT-,  — >  0,  and  t-,/t  — >  oo 

as  a  — >  0.   Since  oo  (p  )  does  not  depend  on  a,  the  validity  of 

each  CD  (p )  can  be  extended  to  cover  the  Interval  i-x^j-x).      The 
n  1 

rigorous  justification  of  this  restriction  will  become  clear 
in  Appendix  VI . 

Substituting  equations  (4.15),  (4.l6),  and  (4.17)  into 
(4.15)  yields  the  following  result: 
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&0 


60 


(?,<)   =  o/^ri.fpj  +y /<7/:fj60ff^^j^j 


(4.18) 


2  ^ 


U- 


:x 


-^ifp^^^fyf[^-f/j^-^si^-'jj-]^[^^o<)j^ 


where  fi  (p)  is  defined  in  (4.15). 

To  utilize  the  concept  of  Hadamard ' s  finite  part  in 
carrying  out  the  limiting  process  a  — >  0  for  equation  (4.18) 
it  is  necessary  to  know  the  behavior  of  cD(r,aj  as  r  — >  oo 
and  that  of  cD(p,a)  as  p  — >0.   It  will  be  assumed  that  cr  cu 
is  the  free-molecule  solution  and  consequently. 


(jo(r  ,a  j 


Jirr 


5 


as 


00,  a 


0 


(4.19) 


And  if  CD(p,a)  and  cD(r,a)  are  to  match  each  other  asymptotically, 
it  is  expected  that 


aj(p,a) 


2a 


5 


^TTP 


3 


as 


0,  a 


0. 


(4.20) 


With  these  results  the  limiting  process  a  — >  0  can  be 
carried  out  for  equation  (4.l8).   It  is  now  clear  from 
equation  (4.l8)  that  a  (a)  =  a    .      Indeed,  divided  by  a    , 
equation  (4.18)  becomes,  in  the  limit  a  — >  0, 


where 


J 
o 


f^o(p) 


^  A  Jx(p) 
3Tr  3  3^^^ 


(4.21) 
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Equation  (4.21)  Is  free  of  a  and  therefore  '^„(p)  is  a 
function  Independent  of  a.  For  this  lowest  outer  solution 
25  (p),  the  boundary  condition  (4.l4)  becomes 


11m  %(p)  =  0  (4.22) 

p  — >00 

It  will  be  shown  In  Chapter  V,  section  V.4  how  the 

existence  of  a  solution  to  equations  (4.21)  and  (4.22), 

2   1 
which  approaches  ^ ^  as  p  — >  0,  can  be  proved. 

With  these  results  for  cd  (p ) ,  the  first  order  outer 

solution  CD-,  (p  )  can  now  be  computed  by  taking  a  co  to  be  the 

free-molecule  solution  and  expanding  Q,{^,   a)    and  K(£-,  — ;  a) 

to  a  higher  order. 


IV.4 .   Inner  Solutions 

In  the  Inner  region  r  <  t,  equation  (4.11)  may  be 
written  as 


T 

^(r,a)  =  fi(r,a)  +  j  K{r  ,n;cL)^{n  A)dn  (^-25) 

1 


.1' 

ax 


where  It  Is  recalled  that 


K(r,  |;  a)^(^.a)d^ 


-  66  - 


dm   f- 


0 

Prom  the  fact  that  r  <  t  in  the  inner  region  and   lim  ax  — >0, 

a  ^  0 
It  Is  seen  that  the  argument  In  the  function  J^  In  the  last 

Integral  Is  very  small.   Asymptotic  expansions  for  the  func- 
tions J  with  a  small  argument  have  been  obtained,  by 
Abramowitz  (1953)-   The  result  for  J^  is 


where  the  asymptotic  symbol  capital  0  has  been  used,  cf. 
Erdelyl  (1956),  p-5)-  Substituting  the  expansions  for  J 
into  n(r,a),  the  following  expression  is  obtained: 


5 


fi(r,a)  =  fi^Cr)  +  afi^(r)  +  ^^^{v)   +  0(a^)n^(r)      (4.24) 


where 


n^(r)  =^rs;.-'^  -]^J^; 


02 


0 
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^^{r) 


•J/n  -pr 


O 


(rcosf  -  Jl-r^Sm^^/sln^j^  c^j 


and   the   function  fi    (r)    Is   the    free-molecule    solution. 

o^  ' 

In  a  similar  manner,  the  kernel  K(r,n;o^)  is  expanded  as 
follows : 


K(r,ii;a)  =  ak^(r,ri)  +  a^K^(r,ii)  +  0{o^)y.^{v  ,r\) 


(4.25) 


where 


Ko(r,ri) 


K. 


<  .  -I  I 


-1  n 
and  b  =  sin   -*■  ,    for   ^  <  r 


^-i 


for   T  >  r 
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For  the  kernel  K(r,  — ;  ex),  again  It  will  be  assumed 


a 

that  r  «  —  and  as  a  result, 
a 


^^^  I,  ^J  =  ^JJ!)(z-i-,[i^,)^0(^')         (*.26) 


After  the  results  in  (4.24),  (4.25),  and  (4.26)  have 
been  obtained,  the  Inner  solution  cd  (r)  and  the  Inner  ordering 
a  (a)  can  then  be  calculated  from  equation  (4.23)  for  n  =  0> 
1,  2,  assuming  that  on  the  R.H.S.  of  equation  (4  23)  to  the 
lowest  order,  cD(ri,a)  can  be  replaced  by  the  free-molecule 
solution.   If  these  are  done,  equation  (4.23)  may  be  put  into 
the  form  of 


''■C 


It  is  recalled  that  as  t]  — >  •x, ,    Q.    {r\)  :^  '^/(tT^^  )    and 
K^(r,ii)*^—  {^   -  ^j—Yl  -  l/r   )  and  on  the  other  hand,  as 
1^0,  oT  (I)  ^  2/371^^.   Consequently, 


T 


lim 
a->  0  ^ 


1 


00 

3  ' 


which  is  a  finite  quantity,  and 

J2[^]^o(4)d^  -  0 
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lim  a 
a->  0    J 


ai 


With  these  results  obtained,  It  Is  clear  from  equation 
(4.27)  that 


C7^(a)  =  a    , 


and 


^ 


CD-, 


-^(r)  =  n^(r)  +   K^(r,Ti)fi^(T])dii  , 


1 


For  the  convenience  of  obtaining  the  second  order  Inner 
solution  cDp(r),  equation  (4.23)  Is  arranged  In  the  following 
form. 


(4.28) 


^7> 


Note  that  (>(J-,  (^1 )  ?^  I/To  as  t)  ->  oo,  and  therefore  the 
integral 


■C 


A 


'^       A  C  ^        / 


converges.   It  should  also  be  noted  that 


DO  OO 


=  O 


Equation  (4.28)  then  yields,  in  the  limit  a  — >  0, 
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^^(r)  =   L(r)  + 


K^(r, 71)03^  (Ti)dri 


A  ,  ,     2 
o^[a)    =  a 


Higher  order  Inner  and  outer  solutions  can  be  obtained 
by  making  alternate  use  of  equations  (^.23)  and  (4.15). 

In  the  next  chapter,  It  Is  shown  that  Inner  and  outer 
solutions  are  Indeed  uniformly  valid  asymptotic  to  cD(r,a) 
and  Go(p  ,a)  . 
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V.   The  Cylinder  Problem--MathematiGal  Foundation 
V • 1 .   Introduction 

The  main  purpose  of  this  chapter  Is  to  prove  the 
Asymptotlclty  Theorem,  I.e.,  to  prove  that  the  Inner  and 
outer  solutions  obtained  In  the  last  chapter  are  Indeed 
uniformly  asymptotic  to  the  true  solution  to  the  cylinder 
problem.   For  this  purpose,  some  preliminary  results  are 
obtained  In  sections  V.2,  V.3,  and  V.4.   After  these  results 
are  obtained,  the  Asymptotlclty  Theorem  can  be  proved  In  V.5. 

The  uniqueness  and  existence  of  the  solution  to  the 
cylinder  problem,  which  are  proved  In  sections  V.2  and  V.3, 
are  not  only  Indispensable  for  the  proof  of  the  Asymptotlclty 
Theorem,  but  are  also  of  Interest  In  their  own  right. 


V.2 .   The  Uniqueness  Theorem 

For  the  rotating  cylinder  In  an  Infinite  gas,  which  is 
at  rest  at  Infinity,  the  mathematical  problem  is  to  find  a 
function  cjD(r,a),  uniformly  bounded  in  1  <  r  <  00,  which 
satisfies  equations  (4.11)  and  (4.12).   The  uniform  bounded- 
ness  condition  is  probably  the  least  restrictive  one  and  yet 
is  sufficient  to  yield  the  desired  results.   For  example,  if 
cD(r,a)  is  uniformly  bounded  for  r  £  [l,°°)  then  co(r,a)  is  con- 
tinuous in  r.   In  fact,  a  stronger  result  has  been  obtained 
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In  the  following  lemma. 

Lemma  5 •  1  •   If  ao(r,a)  is  a  function  uniformly  bounded  In 

1  <  r  <  00  and  satisfies  equation  (4.11),  then  for  r.,  6  [1,N] 

and  Tp  6  [l,N]  there  exist  M(N)  and  &(N)  such  that  I r.  -  r^\ 

I        I  l/2 
<  5(n)  Implies  |cD(r-,,a)  -  cD(rp,a)|  <  M(N)|r^-  rp  |    • 

The  proof  of  this  lemma  Involves  straightforward  but 

lengthy  computations  and  Is  given  In  Appendix  IV.   A  result 

of  this  lemma  Is  that  a3(r,a)  Is  continuous  In  [l,N]  for  any  N. 

Using  this  result  the  uniqueness  theorem  can  be  proved  for  the 

cylinder  problem.   In  the  sequel  frequent  use  will  be  made  of 

the  following  properties  of  the  kernel  K(r,ri;a),  namely 

K(r,T];a)  >  0  (5.1) 

and 

K(r,ii;a)dTi  =  1  -  [fi(r,a)  +  P(r,a;N)]  ,         (5-2) 
J 

1 

where 


^   


and  consequently,  for  any  finite  r  , 
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with  the  result  In  Inequality  (5-3)^  the  uniqueness 
theorem  can  be  proved  as  follows. 

Theorem  5-1  (Uniqueness).   If  v-,(r,a)  and  v^Cr^a)  are  two 
functions  uniformly  bounded  In  1  <  r  <  oo  and  satisfy  equa- 
tions (4.11)  and  (4.12),  thenv-,(r,a)  =  Vp(r,a). 
Proof.   Let  v(r,a)  -   v,  (r,a)  -  v^(r,a)  and  ||v  ||  -      sup  |v(r,a) 
the  least  upper  bound  of  |v(r,a)|  for  r  6.  [1,03).   Suppose  that 
v(r,a)  ^  0.   Then  ||v  ||  >  0.   Equation  (4.12)  requires  that 
|v(r,a)|  — >  0  as  r  — >oo.   Consequently  there  exists  a  finite 
r^  such  that   ||v||  =  Max   |v(r,a)|,  since  v(r,a)  Is  con- 
tlnuous  In  [l,r  ].   And  from  equation  (4.11), 


00 

v(r,a)  =J  K{r,t];a)v{r\,a)(lT]    , 


or 


v(r,a)  I  <  II  V  II 

■J 


r 


K(r,n;a)dT) 


The  last  Inequality  Is  valid  for  all  r  In  [1,°°)  and  in 
particular  for  r  In  [l,r  ].   Taking  the  maximum  value  of 
|v(r,a)|  for  r    [l,r  ]  gives 


or 


V  II  <  k(r^)  ||v 


(1  -  k(r^))  ||v  II  =  0 
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which,  using  (5-3),  Implies  ||v||  =  0  or  v(r,a)  =  0.   This 
completes  the  proof. 


V.3.   The  Existence  Theorem  and  the  Theorem  of 
Continuous  Dependence 

In  order  to  construct  a  solution  to  equations  (^.11) 
and  (4.12),  a  sequence  of  functions  aj(r,a;N)  Is  Introduced 
through  the  following  equations: 

N 
a)(r,a;N)  =  fi(r,a)  +    K(r  ,t]  ;a)cD(ri  ,a  ;N)dT]  (5-4) 

1  I 

Equation  (5-4)  In  fact  defines  a  two-cylinder  problem 
In  which  an  Inner  cylinder  of  radius  1  Is  rotating  at  an 
angular  speed  1  and  an  outer  cylinder  of  radius  N  Is  at  rest. 
This  sequence  a3(r,a;N)  with  Increasing  N  will  be  used  to 
approximate  cD(r,a)  in  the  limit  N  -^   co.   The  physical  meaning 
of  this  approximation  is  self-evident. 

The  uniqueness  of  the  solution  to  equation  (5-4)  for  a 
given  N  can  be  proved  in  a  similar  manner  as  Theorem  5-1- 
The  existence  of  a  solution  cjj(r,a;N)  to  equation  (5-4)  is  a 
consequence  of  the  ,Newmann  series.   In  fact. 


a)(r,a;N)  =  J^ 

n=o 


N 
^   (n) 


^ 


K^'^(r,Ti;a)fi(T],a)dii  +  fi(r,a)      (5-5) 
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is  the  solution  to  i^-^),   where 

N 
K^"^(r,ii;a)  =  J  K(r,t  ;a)K^"~^  ^  (t  ,ti  ;a  )dt  . 
1 

The  convergence  of  the  series  in  (5-5)  follows  from  the 

00 

convergence  of  the  geometrical  series  ^   ["^(N)]  . 

n=o 
A  number  of  the  properties  of  a)(r,a;N)  should  be  noted. 

It  is  clear  from  (5-5)  that  co(r,a;N)  >  0  because  K(r,Ti;a)  >  0 

andn(r>ct)  >  0.   The  fact  that  a)(r,a)  <  1  follows  from  the 

observation  that 

N 
[l-cD(r,a;N)]  -  P(r,a;N)  +    K(r  ,ii  ;a  )  [  l-cD(ri  ,a  ;N)  jdr) , 

^  (5.6) 

where  P(r,a;N)  is  defined  in  (5-2)  and  '|~'(r,a;N)  >  0  for 
^   6  [l^N].   For  then  the  function  [  l-cD(r  ,a;N)  ]  can  be  con- 
structed from  the  Neumann  series  and  is  positive  and  hence 
the  desired  result  follows.   This  result  says  that  a)(r,a;N) 
forms  an  equi-bounded  sequence  of  functions.   Finally,  it 
should  be  noted  that  a3(r,a;N)  is  monotone  in  N;  i.e., 
oo(r,a;N2)  ^   ^^(r ,a;N-j^)  f or  1  <  r  <  N-j^  if  Ng  >  N-^ .   The  proof 
falls  as  a  result  of  the  following  equation  defined  for 
r   6  [1,N], 
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[a)(r,a;N2)-a)(r,a;N-^)  ] 


K(r,r];a)Go(n,a;Np)dTi 


Nn 


N, 


+  j   K(r,ri;a)[co(n,a;N2)-aj(Tl,a;N^)]dT]  , 
1 

which  Implies  that  [a)(r  ,a;Np ) -cD(r  ,a;N-,  )  ]  >  0  for  1  <  r  <  N. 

since  J    K(r,il;a)co(ri,a;N2)clii  >  0. 

^1 
The  following  lemma  Is  an  immediate  consequence  of  the 

fact  that  a3(r,a;N)  is  an  equl-bounded  monotone  sequence. 
Lemma  5-2.   The  limit  function 


00-,.  (r,a)  =   11m  a)(r,a;N) 

N  ->  00 

exists  for  r  6  [ij^^)  and  Is  positive  and  uniformly  bounded 
by  1. 

It  can  be  shown  in  the  following  that  cd-,  .  (r,a)  satisfies 
equation  (4.11). 

Lemma  5.3-   The  function  oo-,  .  (r,a)  satisfies  the  equation 
(4.11);  i.e.. 


co-i  .  (r,a)  -   fi(r,a)  + 
lim^  '  ^     V  '  / 


K(r,T];a)aj,  .  (Ti,a)dTi 


lim 


Proof.   A  limiting  process  N  — >  <»  is  taken  for  equation  (5-4) 


lim  cu(r,a;N) 

N^oo 


N 
==  n(r,a)  +   lim    K(r  ,t]  ;a  )a)(Tl  ,a  ;N)dTi 


N— > 
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and  consequently. 


lim^    ' 


fi(r,a)  +   11m 


K(r,Ti;a)a3  (n,a;N)dTi  , 


where 


CO  (r,a;N)  =  a3(r,a;N)    for   1  <  r  <  N 


=  0 


for   r  >  N 


Since  K(r,ri;a)a)  (r|,a;N)  Is  dominated  by  the  Integrable 
function  K(r,T];a),  It  then  follows  from  Lebesgue  dominated 
convergence  theorem  that 


11m 


r 


K(r,ii;a)a3  (ri,a;N)dTi  =    K(r  ,ii  ;a  )a)^^^(il  ,a  )dTi 


This  concludes  the  proof. 

In  order  to  establish  the  existence  theorem  for  the 
cylinder  problem  It  remains  to  prove  the  following  lemma. 


Lemma  5.4.    lim  cd^  .  (r,a)  =  0. 

—   ^  .     lim^    ^ 

r— >  00 


The  proof  of  this  lemma  Is  a  rather  complicated  affair. 
To  facilitate  the  method  of  proof  It  Is  useful  to  develop  a 
lemma  for  obtaining  more  refined  upper  bounds  for  co(r,a;N). 

For  a  function  f(r),  a  particular  iterate  of  f  is 
defined,  in  b  <  r  <  N,  as 
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b 
iff"]  =  fi(r,a)  +   K(r,Ti;a)a3(T],a;N)dT] 


(5.7) 


N 


+ 


K(r,Ti;a)f  (T])dT]  , 


where  1  <  b  <  N.   This  lemma  for  obtaining  more  refined 
upper  bounds  of  a)(r,a;N)  may  be  stated  as  follows. 

Lemma  5 • 5 •   If,  for  some  e  >  0,  a  function  f (r)  has  the 
properties  that 

(1)   l(f }  <  f(r)   for  b  +  e  <  r  <  N,  and 
(11)   f (r)  >  co(r,a;N)   for  b  <  r  <  b  +  e, 
then  f  (r)  >  ai(r,a;N)  for  b  <  r  <  N. 

Proof.   For  r  <z.    [b,N], 

N 
l[fj-  a3(r,a;N)  =  J  K(r  ,11  ;n)  [  f  (t^  ) -03(1!  ,a  ;N)  ]dii  . 

b 

Hence,  for  r  &    [b+£,N],  It  follows  from  condition  (l)  that 


N 


f  (r)  -  a)(r,a;N)  > 


K(r,ri;a)[f  (r|)-co(n,ci;N)]dri  , 


which  may  be  written  as 


f(r)  -  cu(r,a;N)  >  m 


N 


K(r,Ti;a)dri 


b+e 


b+E 


f 


+ 


K(r,ri;a)[f(n)-a)(r],a;N)]dri 
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where 


m  =  Min   [f(r)  -  a)(r,a;N)] 
b+e,N 


Prom  the  second  property  of  f(r),  It  follows  that 

N 
m  >  m 

b+e 
which  implies  ni(l-"p(N))  >  0  or  m  >  0,  and  consequently. 


r 

K(r,ii;a)dTi  , 


f(r)  -  a)(r,a;N)  >  0   for   r  6  [b,N]. 

This  completes  the  proof. 

As  a  special  case  for  b  =  1,  e  =  0,  Lemma  5-5  reduces 
to  a  theorem  of  Snyder,  cf.  De  Marcus  (1957)- 

Using  Lemma  5-5  it  can  be  shown  in  Appendix  V  that 

f(r)  =  c/r  is  an  upper  bound  of  cD(r,a;N)  for  all  N  and  a, 

where  c  is  a  constant  independent  of  N.   It  then  implies 

that   lim  co^  .  (r,a)  =  0,  which  is  Lemma  5-^- 
.     lim^  ^    I  ' 

r  ->  00 

The  results  of  Lemmas  5-3  and  5-''+  can  be  summarized 
in  the  following  existence  theorem. 

Theorem  5-2  (Existence).   There  exists  a  solution  co-,  .  (r,a) 
i' ^  lim^    ' 

to  equations  (4.11)  and  (-'1.12). 

It  is  seen  from  Theorem  5-1  that  oo,  .  (r,a)  =  a)(r,a). 

•^         lim^  '    '  \    '    I 

After  the  uniqueness  and  existence  theorems  are  estab- 
lished, it  becomes  a  simple  matter  to  show  that  the  solution 
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a3(r,a)  to  the  cylinder  problem  depends  continuously  on  the 
boundary  term  fi(r,a).   A  more  precise  statement  can  be  given 
In  the  following. 

Theorem  5.3  (Continuous  Dependence).   If  the  function  g(r,a) 
satisfies,  for  0  <  e  <  1, 


:(r,a)  =  G(r,a)  + 


K(r,T];a)g(Ti,a)dTi  , 


where 


G(r,a)  -  fi(r,a)|  <  sfi(r,a). 


then 


g(r,a)  -  a3(r,a)|  <  eco(r,a) 


Proof.   Let  the  functions  e(r,a)  and  E(r,a)  be  defined  as 


e(r,a)  =  g(r,a)  -  cD(r,a) 


E(r,a)  =  G(r,a)  -  fi(r,a) 


Then,  e(r,a)  satisfies  the  following  equation. 


e(r,a)  =  E(r,a)  +    K(r /n  ;a  )e  (rj  ,a  )dT] . 

1 


And,  according  to  theorems  5.I  and  5.2,  there  exists  a  unique 
solution  to  (5-8)  defined  by 
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Therefore , 

N 

Co    r      (nJ 


which  completes  the  proof. 

In  passing.  It  may  be  pointed  out  that  all  the  results 
obtained  In  this  and  the  last  sections  are  valid  for  any 
arbitrary  values  of  a. 


V.^.   The  Outer  Solutions 

In  the  section  to  follow  the  powerful  Theorem  5-3 

will  be  used  to  prove  that  the  Inner  and  outer  solutions 

'Zoo:)     and  Z  a  cci  are  Indeed  uniformly  asymptotic  expan- 
nn        nn  j        j    f  t- 

slons  for  the  solution  co(r,a)  for  a  «  1.   In  order  to  do 
this,  it  is  necessary  to  show  that  there  exist  functions 
co'  satisfying  the  n-th  order  outer  equations.   If  cd  exists, 
the  CD  can  be  obtained  by  quadratures.   It  is  also  to  be 
shown  that  cu  (p)  have  the  expected  asymptotic  behavior  as 
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P  ->o. 

One  way  to  show  that  co  { p )  exists  Is  to  use  the  Fourier 
transform  technique,  cf.  Grad  (1959)  and  Sirovlch  (i960), 
and  prove  that  the  outer  equations  are  equivalent  to  the 
Fourier-transformed  equations  with  proper  source  terms.   Then 
the  existence  of  the  solutions  to  the  Fourier-transformed 
algebraic  equations  Implies  the  existence  of  cd  (p )  .   The 
behavior  of  co  (p )  for  p  — >  0  can  then  be  obtained  by 
explicitly  evaluating  the  solutions  obtained  by  the  Fourier 
technique  for  p  — >  0. 

It  Is,  of  course,  highly  desirable  to  secure  a  direct 
proof  for  the  existence  of  a  solution  to  the  outer  equations. 
Although  such  a  direct  proof  has  not  been  found.  It  will  be 
Indicated  In  what  follows  a  possible  way  to  arrive  at  this 
proof  and  the  difficulties  Involved  therewith.   For  example, 
consider  equation  (4.21)  for  oj  (p).   If  a  function  h  (p)  Is 
Introduced  through  the  equation 


where 


K^°^  {p  A)    =   K(p,^),  and 


K^"^P,^) 


00 


K(p,t)K^"-^)  (r,4)dt  , 


o 


then  equation  (4.21)  may  be  written  as 
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h,(p) 


H^(P)  + 


K(p,^)h^(^)d^ 


o 


where 


H,(P) 


^  K^^hp  ,^)'^A^)d^    , 


^o 


(5-10) 


which  Is  uniformly  bounded  for  0  <  p  <  oo. 

The  boundary  condition  (4.22)  goes  over  to 


11m  h  (p)  =  0 

p  ^00 


(5.11) 


as  a  consequence  of  the  fact  that  as  p  — >  0, 


(^ 


1    <^ 
n^o  J^ 


K(")(p,^)n^(^)d^]— >0 


Equation  (5. 10)  Is  different  from  equation  (4.21)  In 
that  the  Inhomogeneous  term  H  (p )  is  uniformly  bounded 
while  Q,    (p)  diverges  as  p  — >  0. 

It  can  be  shown  In  a  manner  similar  to  Theorem  5-1 
that  If  there  exists  a  uniformly  bounded  solution  to 
equations  (5-10)  and  (5. 11),  then  the  solution  Is  unique. 
Furthermore,  It  can  be  shown  by  constructing  the  Neumann 
series  that  there  exists  a  uniformly  bounded  solution 
^qCPjN)  to  the  following  equation. 


hQ(p,N)  =  H^(p)  + 


N 
K(p,e)h^(^,N)d4 


(5.12) 
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If  a  function  T(p)  can  be  found  such  that  for  0  <  p  <  ~, 

00 

Ho(p)<T(p)  -  [  K(p,nT(^)d4  ,  (5.15) 


then,  h  (p,N)  Is  equl-bounded  and  consequently  the  existence 
of  a  uniformly  bounded  h  (p )  to  equations  (5-10)  and  (5-11) 
will  follow  in  a  manner  similar  to  the  proof  of  Theorem  5-2. 
Then  cu  (p  )  has  the  expected  behavior  as  p  — >  0,  as  can  be 
seen  from  the  definition  of  h  (p ) .   However,  It  has  not 
possible  to  find  the  function  T(p)  satisfying  (5.13)- 


V.5.   Uniformly  Valid  Asymptotlclty 

In  this  section  the  powerful  theorem  5-3  Is  used  to 
prove  that  the  Inner  and  outer  solutions  obtained  In  the 
last  chapter  are  Indeed  uniformly  asymptotic  expansions 
for  the  solution  cD(r,a)  for  a  «  1.   It  will  be  assumed 
that  the  existence  of  co  (p),  which  have  the  expected  asymp- 
totic behavior  as  p  — >  0,  have  been  proved  In  a  manner 
indicated  in  section  V-^- 

Consider,  for  example,  the  function  g(r,a)  defined  by 


y^*•    .    .-<%-* 


g(r,a)  =  a^(a)oj^(ar)     r  >  t  , 


Oo(cx)a)^(r)     r  <  T  . 
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and  define 

e(r,a)  =   g(r,a)  -  co(r,a) 

E(r,a)  =  G(r,a)  -  n(r,a)  , 
where  „ 

00 

/ 

G(r,a)  =  g(r,a)  -   K{v ,r];a)g{n,a)dn    . 

J 
1 

In  order  to  prove  that  for  some  M, 

|e(r,a)|  <  /aMa)(r,a)  ,    as  a  -^  0, 
It  Is  only  necessary  to  show  that 

|E(r,a)|<yaMfi(r,a),    as  a— >0,        (5-1^) 

for  then  the  desired  conclusion  follows  linmedlately  from 
Theorem  5-3- 

The  proof  for  inequality  (5-14)  Is  very  lengthy  and 
may  be  found  In  Appendix  VI.   Similar  results  may  be 
derived  for  higher  order  asymptotic  solutions  by  following 
the  method  used  In  Appendix  VI,  although  the  proof  will 
become  increasingly  tedious. 
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VI.   Discussion 

VI. 1.   The  Asymptotic  Theory 

A  uniformly  valid  asymptotic  theory  of  nearly  free- 
molecule  flow  has  been  developed  In  Chapter  III  for  the 
linearized  Boltzmann  equation  for  hard-sphere  or  power-law 
molecules.   In  the  last  two  chapters  the  theory  has  been 
successfully  applied  to  the  simple  cylinder  problem  based 
on  the  BGKW  equation  and  the  approximate  solutions  so 
obtained  have  been  proved  to  be  asymptotic  to  the  true 
solution  throughout  the  physical  space.   The  success  of  such 
a  proof  for  this  specific  simple  problem  may  be  taken  as  an 
Indication  of  the  correctness  of  the  general  asymptotic 
theory,  although  an  Immediate  extension  of  the  method  of 
proof  cannot  be  found  to  cover  the  general  nearly  free- 
molecule  flow  problems. 


VI. 2.   Drag 

For  a  uniform  flow  past  a  two-dimensional  cylinder 
of  arbitrary,  convex  cross-section,  the  aerodynamic  drag  D 
Is  found  to  be  (cf.  section  III. 5), 


2 
D  =  D^  j^  +  a^n(l/a)D-j^  +  aDg  +  0(a  in  a)  .     (6.1) 
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It  should  be  noted  that  the  first  order  correction  D^ 
Is  due  to  the  contribution  of  colllslonal  effects  In  the 
tall  of  the  inner  region  and  the  beginning  part  of  the  outer 
region.   In  other  words,  D-,  is  caused  by  the  collisional 
effects  in  the  region  where  inner  and  outer  solutions  over- 
lap.  In  this  overlap  region  the  free -molecule  solution  is 
still  the  lowest  order  description  of  the  flow  field.   This 
may  also  be  seen  from  the  expression  for  D-,  ,  in  which  the 

function  G   is  related  to  the  free -molecule  solution  5'^ 
o  o 

through  the  limiting  process 

G  (c,  e)   =  iim  [^  ^,{^i.^,0;r^B^^)} 

It  is  then  concluded  that  D-,  can  be  calculated  by 

Willis'  method  and  for  this  purpose  only  the  knowledge  of 

vp  is  needed, 
-^o 

^  ^^ 

Both  $  and  ^  are  needed  in  order  to  compute  the 
•^  o     -"-o 

next  order  correction  term  D^ •   The  inner  solution  $  is 
simply  the  free-molecule  solution  "f  .   For  the  outer 
solution  ^  it  is  necessary  to  solve  equations  (5.21)  and 
(3.22)  numerically.   While  it  is  true  that  equations  (5-21) 
and  (5.22)  may  not  be  much  easier  to  solve  numerically 
than  the  original  problem.,  i.e.,  equations  (2.l6)  and  (2.17), 
the  advantage  of  the  asymptotic  theory  lies  in  the  fact 
that  equations  (3-21)  and  (3-22),  unlike  equations  (2.l6) 
and  (2.17),  are  independent  of  the  parameter  a  and  can 
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therefore  be  solved  once  for  all  small  a 


¥1.3-   The  Lowest  Outer  Ordering 

One  of  the  most  Important  quantities  in  the  theory  Is 
the  lowest  outer  ordering  a  (a)  which  determines  the  order 
of  magnitude  of  the  outer  region  effects  on  the  aerodynamic 
forces.   It  appears  therefore  worthwhile  to  examine  the 
physical  origin  of  this  outer  ordering  without  going  through 
the  whole  mathematical  structure  of  the  theory.   Consider  a 
two-dimensional  problem  in  which  the  mean  free  path  l/a  is 
used  as  the  length  scale .   To  the  outer  solution  ^,  the 
contribution  due  to  the  inner  solution  can  be  neglected  in 
comparison  with  that  due  to  molecules  coming  directly  from 
the  boundary  of  the  body.   Therefore,  to  the  lowest  order, 
the  outer  field  Is  created  entirely  by  the  boundary  contribu- 
tion as  viewed  from  a  distance  much  larger  than  the  body 
dimension  away  from  the  body,  i.e.,  p  >  aT.   Then  in  the 
limit  a  — >  0,  the  outer  flow  field  is  created  by  a  source 
at  the  origin  which  now  is  the  only  quantity  in  the  problem 
which  involves  a. 

Because  of  the  linearity  of  the  problem,  the  magnitude 
of  the  outer  field  depends  solely  on  the  strength  of  the 
source.   The  lowest  order  expression  for  the  source  term  may 
be  obtained  for  hard-sphere  molecules  as  follows. 
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If! 


Co   go  "® 


which  determines  the  outer  ordering  cr  (a)  =  a. 

Note  that  the  ct/p  dependence  Is  of  a  purely  geometrical 
nature,  resulting  from  the  <J5  Integrating  In  (6.2)  where  the 
angles  (H)-.    and  (ff)     subtended  by  the  body  at  {p  ,6)   are  both 
of  the  order  of  a/p . 


VI.^.   Three-dimensional  Problems 


In  view  of  the  geometrical  nature  of  the  lowest  outer 
ordering  a  (a)  the  order  of  magnitude  of  the  outer  solution 
may  be  estimated  for  three -dlm.enslonal  problems  without 
going  through  a  detailed  mathematical  analysis.   Since  the 

angles  subtended  by  a  three-dimensional  body  Is  of  the 

o   o  ,  p 

order  a  /p  ,  the  lowest  outer  ordering  o\,(ct)  will  be  a  . 

Prom  this,  together  with  the  fact  that  the  free  molecule 

solution  decays  as  l/r  at  large  r.  It  may  be  concluded 

that  the  drag  In  a  three-dimensional  problem  has  the 

expression 


2 

D  =  D^  +  aD-j^  +  a  in(l/a)D2  +  a  D^ 
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where  D-,  and  Dp  are  entirely  due  to  the  collislonal  effects 
in  the  inner  region  and  may  be  computed  by  Willis'  itera- 
tion method.   For  the  next  order  correction  D^,  again  ^^ 

I)  -^o 

must  be  obtained  by  solving  the  corresponding  equations 
numerically. 


VI. 5.   Non-Linear  Problems 

Although  the  analysis  in  the  present  study  is  restricted 
to  linear  problems,  some  of  the  results  may  be  extended  to 
non-linear  cases,  provided  the  body  is  not  cold  and 
M/k  «  1.   For  example,  consider  a  two-dimensional  non- 
linear problem.   The  flow  field  outside  of  a  fictitious 
circle  C  of  radius  x  will  be  of  a  linear  nature  provided 
T   is  taken  to  be  very  large  compared  with  the  body  dimen- 
sion.  Then,  the  asymptotic  theory  is  applicable  to  the  flow 
problem  outside  of  the  fictitious  circle  C  if  the  correct 
boundary  value  can  be  prescribed  on  C.   The  radius  t   should 
be  so  taken  that  "^ J '^   — >0,  t  ^^ooasa->0.   Again  the 

lowest  order  outer  solution  ^  may  be  obtained  for  p  >  ar, 

«^ 
if  the  source  term  '^  Is  calculated  from  the  correct  flcti- 

tious  boundary  condition  on  C.   But,  since  V'^p,  ~^   °°  ^-^ 

a  — >  0,  then  to  the  lowest  order  the  function  ^  can  be 

calculated  from  the  real  free -molecule  boundary  condition  on 
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the  body.   In  this  manner,  ^  can  be  determined  and  has 
precisely  the  same  form  as  In  the  linear  problems  and  It 
Is  clear  that  o'  (a)  =  a.   Now  the  Inner  solutions  cannot 
be  obtained  by  the  linear  theory.   However,  in  the  inner 
region  for  1  <  r  <  t,  the  colllsional  effects  can  be  con- 
sidered as  a  perturbation  and  Willis'  method  applies.   It 
must  be  noted  that  the  former  expression  for  the  scattering 
probability  distribution  P(c_,c^;  D(r,0))  (see  section  II. 3) 
is  not  valid  for  non-linear  problems.   The  usual  non-linear 
free -molecule  boundary  condition  should  suffice  for  the 
purpose  of  finding  the  lowest  order  outer  solution  ^^. 
Again,  the  drag  has  the  expression 

D  =  D^j^  +  ain(l/a)D^  +  a^^ 

and  D  can  be  calculated  by  quadratures  using  the  tail  of 
the  free-molecule  solution.   To  the  next  order  Dg  may  be 
computed  from  $  and  the  distribution  function  f  in  the 
inner  region  obtained  by  Willis'  method.   Analogous 
results  may  be  obtained  for  three-dimensional  problems. 
In  concluding,  it  must  be  noted  that  the  condition 
M/k  «  1  may  be  too  restrictive  for  high-speed  aerodynamic 
problems.   A  complete  non-linear  analysis  is  necessary  to 
remove  this  restriction. 
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Appendix  I.   Diffuse  Boundary  Condition 

The  scattering  probability  P(c_,c_-,;  D(r,0))  Is  derived 
here  for  the  diffuse  re-emlsslon  boundary  condition. 
Let  the  re-emlsslon  distribution  function  be 


T. 


3/2    (  ± 

\=   ^h^2W^^         ^M-    ^^-T^  hn>  ° 


Then,  the  conservation  laws  can  be  expressed  as 


///  r/.  ^^c  =  fff  /-^  n^f  ;^nc/3c 


:i-i) 


On  carrying  out  some  of  the  integrals  and  neglecting 

2 
terms  of  U  and  the  product  terms  of  U  and  f ,  the  above 

equations  become 


Let  n,  and  T,  be  linearized  according  to  the  following 
formulae : 
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n,  =  n  ( 1  +  n,  ) 
b    o^     b' 


^b  =  ^o^l  +  ^b) 


Then  equations  (l.2)  become 


r^^  ^^r,'-  - Jf  Li C..9.  -f/[C/'"^  ^^'^ 


and 


^  3_r.^=-^LJcas^n^n^Jn^  c.g.^^^c 


^b  -^iTb 


.-^  .^r.^ 


C^<6 


I        I 


from  which  n,  and  T,  can  be  determined;  In  fact 


n,'=  -  Iff  c7  co.e,  +  ^  \{\i'U  -  ?>„  f  d  ^c  ^ 


c^^o 


i-  ^  2 


V  =  -fuc..e.t4///eV|-'JC.fd'c. 


C^<o 


Henceforth, 


^  -    ^c    ^2rrfiTo 


(: 


I 


(1-3) 


/-^  n. 


2 


-fr/^c^rj^ 


-c 


I        \3/: 
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^  =   %   +    (c^-  |)T^    -   2c-U 


TT  /.  .   C 


2C.U    -   ^J!(,+  ^yL)c,se 


0 


PC^.^^-D^,;;  =  A  e-\34^'-c^c^>„ 
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Appendix  II 

Now  the  asymptotic  expansion  (3-17)  will  be  derived 
for  ^  (c_;  ^,  0;  a)  in  the  outer  region. 


a 


i^fi,C.Or..,)^-A    /c..^      nl  (Il-l) 


'k^'-'L^.^}='j^[i(ijf,»;^)^i('-jtK^)] 


where 


--^,(W  «^   c 

and  K-,  (c_,c_-,  )  and  Kp(c_,c_,  )  are  defined  in  (2.5). 

A  number  of  ordering  functions  A-,  (a )  and  Ap  (a )  are 
chosen  with  the  properties  that  as  a  — >  0  :  A-,  — >  0, 
Ap  — >  00,  VAp  — >  0°  and  A-.T  — >  CO.   Then  for  s  <  Ap , 
1//?  <  s-i  <  T,  c  <  A„,  c„  <  t  and  |s-s-,  |  >  A,  ,  the  kernel 
Kp(c_,c_,  )  can  be  expanded  as 


and 


exp{-|;(fc«.f.-c^fc^^;}  ^   ex|:>f-il(S|z£2_p| 
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Therefore  the  f  unction  "^P  (c_;p  ,0  ;a)  may  be  written  as 
\{'^j?.^j'^)  ^    I  T-^.T:;  i:\   S-A,;K^J.,c,;ex^C_|;^)]   (II. 2) 

In  the  limit  a  — >  0,  the  last  three  integrals  of 

(II -2)  all  vanish  and  the  sum  of  the  first  two  Integrals 

-v/s^  p 
converges  to  l[oo,-oo;  0,<»;  Kp-^(c_,c_  )e       ].   As  a  result, 

it  is  obtained  that 

-I2(c_,p  ,e;a)  ^-^2(c_,p  ,0;C)  . 
Similarly,  as  a  — ^O, 

-fe^(c,p,0;a)  ;^-R^(c,p,0;O)  . 
Summarizing  the  results  gives 


f^'-.l.».^)  ^  -^[4,(^-.ffi;o)^l(s,f.ejOj] 
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which  gives  the  explicit  expression  for  (3. 17)  as 

$q(c;P,0)  =  2[^^(c;p,0;O)  +  i^{^;p  ,Q ;0)] 

The  derivations  in  this  appendix  can  be  made  rigorous 
by  using  inequalities. 
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Appendix  III 

The  derivation  of  expression  (4.15)^ 

p       r   f 

0 

where  a/p  <  1/t. 


A  function  T-,(a)  Is  Introduced  with  the  property  that 

2 
11m  T-,  =  0,   11m  t/t  —>o:>,      and   11m  a  /t  t  =  0.   The 
a^o         a^o  a->o 

integration  in  J^  will  be  split  into  two  ranges. 


■+  }s^-^f  j-S^-i(i^Ji^-!^^^fs^f)\^S 


T^i 


The  first  Integral  in  the  last  equation  is  always 

exponentially  small  because  the  condition   lim  t/t-,  =  oo. 

a  ^o 

And,  in  the  second  Integral  of  the  last  equation  the 
exponential  function  can  be  Taylor-expanded.   The  result  is 
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and  consequently  In  the  limit  a   — >  0, 

Higher  order  expansions  can  be  obtained  without 
difficulties  by  taking  more  terms  in  the  Taylor  expansion 
for  (III. 2) . 
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Appendix  IV 

In  this  appendix,  a  proof  Is  given  for  Lemma  5-1- 
Lemma  5  •  1  •   If  a)(r,a)  Is  a  function  uniformly  bounded 
In  1  <  r  <  00  and  satisfies  equation  (4.11),  then  for 
r-j_  e  [1,N],  r^  £  [1,N],  and  1  <  N  <  ~,  there  exist  M(n) 
and  6(N)  such  that 


I     I      I  1 /2 
\ud{r^,a)    -   cD(r^,a)|  <  Mlr^-r-j^l  / 


Proof.   Assume  rp  >  r, ,  and  write 


cD(rp,a)  -  a)(r  ,a)  =  n(rp,a)  -  n(r-,,a)  +  I 


where  I  =  I-,  +  Ip  +  I^  +  I^, 


^1  =A,[^-'i|;r,s;.e,r,^.Ojr,};  -] 
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and 


The  following  results  are  obtained. 


(a)      ll^l    <    \l^^\    +  1^2  +  I13  +  Ii4 


where 


Hill  =|.i(-'JAji^j  s- 1^- ^  J  '"'r^  -J 


•12 


■15 


■14 


and 
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From  the  fact  that  co(t])  <  1 ,  It  can  be  shown  that 


^iil  <  ^y.[o)lr,''-ff'-^2J^(o)r^il^--r;')(Lj-^^) 


^    ^JJ^)  r^^\K-r, 


% 


For  I-ip'  ^1^'    ^"^  ^ik'    ^^    ^^   found  that 


Il2  <  J2(0)  ^^   (^2-^1)  <  2J2(0)(r2-r,) 


-13  <  (r2-r^)j2(0)(2r2+  l) 


and  finally. 


^14  <  IIfo)[(lK,tiXrl-vJ+TT  (r,^-ri^J%  ^' 


^    (2r;  -r  ^Ty)X(c.)  |r,-in 


Jl 


Hence , 


1^1    <   C-^r^   J2(0)|r2  -r^|^   ^ 


103   - 


whe 


re  the  constant  C-,  Is  Independent  of  r-,  and  r^ 


(b)   ll^l  <  I21 


+  I 


22  ' 


where 


-21 


and 


'22 


=  A,[».i  J  r,/.;  I   -] 


It  is  found  that 


-h 


^r^O 


and. 


I22  <  J2  ^  ° )  ^  ^2  "  ^1  ^  ^  ^  +  1^2  ^  -^^  ^1  ^ 


And,  therefore. 


I2  I  <  M2  I  rg  -  r^ 


1/2 


where  M2  is  a  function  of  N  and  a 
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(c)   I  1^1  <  I^^  +  I32  +  I33  ' 


where 


I3I  =  I  Z  (-O7V,  [s-^^fj^S-^e.OOj  «f^)  ;  +J 


I52  ^A2t"'t.5-"frj^^'"e.^;  ij  +  1 


and 


^33  =A|L^''"''l^,'^;^^-^^.^-^j  L-h] 


The  results  are: 


I31I  <  3^1  X  i^)W.-rS\  ^J.(o){f-,-Y,) 


< 


^r.^ 


[zr,.  jf^    JlC-r, 


»/. 


and. 


J  (0)  (r  -r  ) 

T     /  _ii ±— 

^32  ^   a    r^r^ 


\r2        2' 
yrp  -  r-i 

I35  <  -^H?7^  Jjt") 


Hence 


-3 


<[..,.  IJ^jK-, 


/^ 
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(d)     |I^|    <   1^^  +   1^2  +   I43    ' 


where 


•41 


u"=l  ^ 


and, 


•^ 


-43 


=  A 


^^^^)°^'^'\;   ': 


-] 


The   results   are    found   to  be: 


hi< 


[.r^^,C.ir.)J.(^)^^^L(^j\ 


r^■r^ 


A 


I42  < 


lor^]  Ir^-r^l 


y^ 


and 


l43< 


[i  ^/^^^J^rr.-ojj^roji^-^l 


'/z 


Hence 


II4I    <  M^(N)|r2  -rJ^/2 
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where  Mi^(N)  is  a  function  of  N. 

Finally,  It  can  easily  be  shown  that 


^2 


Collecting  the  results  obtained  In  the  above  gives 
|a3(r2,a)  -cu(r^,a)|  <  M(N,a  )  |  r^  -  r^  |  ^/^  , 
which  proves  Lemma  5.1. 
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Appendix  V 

It  will  now  be  proved  that  f(r)  =  c/r  for  some  constant 
C  Is  an  upper  bound  for  cu(r,a;N)  for  all  N. 
Proof.   Consider  the  difference 

b 
l[f j  -  f  =  fi(r,a)  +  I  K(r^,n;ct)a3(T],a;N)dTl 

1 

-  f(b)n2(r,a)  -  f(N)fi^(r,a) 


where 


+  T^  +  Tg  +  T^  +  Ti^ 


i/V^ 


^2 


(r,a)  =^j  jj',[o/(rCos$-j7^T^si^)]^^^^ec^e 
i   c 

4  /  Sjo((rc^s6-^Jb'^-r'-s!^^e}1^'^''9^ie  j 


o   ^ 
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T 


3 


=  T[krj   i/'n-'-f  ^Sm  '^  J   i'irj)    .   -f-  ] 


T,  v'T'ft.-^^^/n-i^.^;^),-] 


and. 


i     d 


7[f->;r,cij3r^);  tj  =J=/'^p&5r^JS/n''^[o^|rc«7<S^iy^^-r'j;n^^|] 


Ci.   c 


If  the  constant  c  Is  such  that  f(b)  >  1,  then  equation 
(V.l)  may  be  written  as  an  Inequality 


^f]  -  f  <    (K,r,r|;a)aj(T),a;N)dTi  -  f(N)fi^(r,a) 


that 


and 


+  T^  +  T^  +  T^  +  T^  . 


It  will  now  be  shown  that  there  exist  b  and  ^  such 


l(f]-  f  <  0    for    b+e  <  r  <  N  (V.2) 


f  >  co(r,a;N)    for   t)  <  r  <  r+e 


Then,  It  follows  from  Lemma  5-5  that 


a)(r,a;N)  <  f(r)  =  c/r  . 


In  order  to  prove  (V.2),  the  functions  T. 's,  1  =  1,2,3,4, 
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will  be  evaluated  for  r  >  b+e,  where  b  »  l/a.   To  this 
end.  It  Is  noted  that 

and  accordingly,  T.  <  0  and  T^  <  0  . 

Now  consider  the  function  Tp •   With  a  suitable  change 
of  variables,  Tp  may  be  written  as 


o     on     qo  ^ 


where 


Tot  =-§(1  P'^^^^PM'T^JTr -^^-^cdSc^Q 


It  should  be  noted  that  the  Integral  of  the  function 
Tp.  contains  r  only  in  the  exponential  function.   There- 
fore, when  ar  is  very  large,  the  asymptotic  value  of  this 
integral  can  be  evaluated  by  the  Laplace  method,  cf . 
Erdelyl  (1956) .   The  result  is  found  to  be 
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For  Tpp,  it  can  easily  be  shown  that 


C      O 


<   f(^jQj(r,^J 


If  e  Is  taken  to  be  equal  to  b,  then  the  Leplace 
method  can  be  used  to  evaluate  T^,  .   The  result  Is 


■h  -   ^'-  I  S  ^r^        IG   ^Y-^ 


> 


And  for  r  >  2b,  it  can  be  shown  that 

^  r 

5  0  U  [f^) 

This  result  Illustrates  the  necessity  of  Introducing  £ 
into  Lemma  5-5- 

Summarizing  all  the  results  obtained  above  implies 
that  there  exist  b  for  which  ab  »  1  and  e  =  b,  such  that 

i{f  I  -  f  <  -  ffM^o(,-~.j  ^  ouf^]) 


-  Ill  - 


Henceforth, 


If   -  f  <  0   for   r  >  b+e  =  2b 


The  constant  c  Is  taken  to  be  c  =  2b;  then 

f (r)  =  c/r  >  1  ,  b  <  ^  <  b+e  =  2b 

or 

f (r)  >  a)(r,a;N)     for    b  <  r  <  b+e  =  2b. 

Then,  according  to  Lemma  5-5j 

Pb 

-^  >  a)(r,a;N)    for   b  <  r  <  N. 

The  analysis  here  is  valid  for  any  N,  and  hence 


Ph 

Y-   >  CD(r,a;N)    for    all  N, 

or, 

2b 
r 


>  00-,  .  (r,a)  , 

which  implies  ao-,  .   (r,a)  — >  0  as  r 
^      lim  ^    ^   ^ 
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Appendix  VI 

In  this  appendix  the  inequality  (5.l4)  will  be  proved; 
that  is,  for  some  constant  M, 

|E(r,a)|<yaMfi(r,a)    as    a^O. 

This  will  be  proved  for  the  outer  region  and  the 

inner  region  separately.   The  function  T(a)  is  taken  in 

-1/2 
this  proof  to  be  a  ' 

(i)   In  the  outer  region,  p  >  qt,  and 

where 

CO 


II  =  o^'i[l<(^J;U)-^KC^!)]:0,(^)cil 


2  =  I  ^(Ip^)^o(V'^) 


I 

I 


h  ^  "^^j  i<^f4)^off)^l 


c 


1^  -o/^|/<r/>.j;^orjj^5 
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and 


-J  "C, 

T-,  =  a  /  . 


From  the  definitions  of  K(r,n;ci)  and  K(p,|)  it 
follows  that 


and  consequently 


r^i 


=  0('^Vs}[q(f.s)  +  T.(f)jr^._ 


< 


0(^^) 


f 


T    Z^ 


Therefore , 


I,<  /^  M.a^n^(p) 


1^  "  o' 


(VI  1) 


For  the  Integrals  1^   and  I^,  It  is  noted  that  I  «  p 
and  rj  «  p/a ,  and  consequently 


/rr^^^;^;^cr^V-^v^;^ 


^ 

J 
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^-v^ 


and 


K(f,j)  ~  0(0  ^r ' 


Then,  if  it  is  recalled  that 


A 

GO 


,3 


O 


03 


(t)   <   c-,/11^  for   some   c 


o 


(^)   <    cV^^ 


CO 


1  I 

for   some   Cp    ,  j  -»  o 


and 


Jo(p)     ^^ 


the  following  results  are  obtained  for  Ip  and  I^,  namely, 

dT)  (VI.  2) 

a^Jo(p) 


P    1 


<  O(aT^)  j- 


<  /a  M^  a^fi^(p) 


I^  = 


3  JpCp)  ^ 

P    ^ 
o 


de 


a^Jp(p) 


=  0(aT  )  y 

P 


<  /a  M^a^fi^(p) 


(VI. 5) 
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For  the  Integrals  Ij^  and  I;-,  the  functions  cd^(^), 
K(p,l),  and  k(^,  n^  ct)  can  be  simplified  by  noting  that 
^  «  1  and  ari  «  1,  and  the  function  o)  (ii)  can  be 


simplified  for  t]  »  l .   The  results  are. 


.. ;  M:  ^^(M^dy^  '^'^ 


-  J 


4/ 


Of<^')    Z 


^Z 


For  Ihp  the  order  of  Integrations  can  be  inter- 
changed to  obtain  the  following  results: 


-^^2. 


J--VC 


'-'/^z[f''i][^s/.p 


I 


-  ocj 


jOS/h^ 


^/'M-^r,_y;77v^ 


/ 


r 
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which,  for  p  <  a  '  ,  reduces  to 


,?,         ^  Jp(p)P     ^ 

P 


;5  J2(P 


3 


and  for  p  >  a 


Jo(p) 


2_2 


Or^^Jl^2<0(a5)^[l-(l-l«^- 


2    2 
P 


)] 


<  0(a5)  :^a2T2 

P 

<  /^  M|a^fi^(p)  . 


For  Ir-j  the  results  are 
D 

*-'    J  ^v^.  -;r"  3'*^  \ 


JT^r^^// 


^ 


-f-  o  r^ ;  /    + 


x,  ^  ^  J)' 


7  J^7-^^/.> 


^  X^/    Y-  or^j  Z^:^ 


'4  J/^5^  ^        I 


and 
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S^-'^  , : ,   .   ,.^  i 


I 


■s^ 


-  omjF^}li4/.t[J^-'-^_ 


^Zsm.  ^j^ii^ . 


r 

By  separate  consideration  for  p  <  qt^  and  p  >  ar^,    it 
can  be  shown  that 


3; 


0{oi)l^2  <   ^"^  M,a-^fi^(p) 


Finally,  note  that 


^Z  ^In-^ 


^TM^o^'^i^o^f ) 


Therefore, 


|I^  +  I3I  <  Ka  M,a\(p)  (VI.^) 

Summarizing  (VI. l),  (VI. 2),  (VI. 3),  and  (VI. 4)  gives 
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H(£ocj|  <  rT^M,  ^^QoCf) 


from  which  it  follows  that 


^C^^^)\^!^M-^(loiJ 


(VI. 5) 


(ii)   In  the  inner  region,  r  <  t. 


E(r,a)  =  oj  (r)  -  fi(r,a) 


K(r,ri;a)g(ii,a)dTi 


a)^(r)  -  fi(r,a)  -  (ig  +  I^) 


where 


,T 


-6 


K(r,Ti;a)a^^(T])dTi  , 


and 


Iry  =  a 
'      J 


K(r,  ^;  a)a3  (4)d^  . 
a     o   ^ 


1/2 


Since  r  <  a  ^  ,  it  is  seen  that 


Ig  =  a   K^(r,T])<i^(T])dTi  +  0(a^: 


Ki(r,T])^^(r|)dTi 


or 


Igl  <  a    K^(r,T])a3^(n)dii  +  0(a) 
1 


K-^(r,ii)a3^(ri)dri 
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On  account  of  the  facts  that  as  r  — >  oo. 


^^(r) 


37rr 


3 


and 


KQ(r,Ti)^^(n)dT] 


K^(r,Ti)2)^(Ti)dri 


2 


it  may  be  concluded  that 


Igl  <  /a  M^fi^(r) 


(VI.  6) 


And  I„  may  be  written  as 


ffo 


=  X„  -t  T 


72 


where  t^  Is  taken  to  be  a 


-V5 


Then  for  !„-,  the  kernel 


K(r,  — ;  a)  and  co  (2,)  may  be  simplified  by  noting  ar  «  1, 
^  «  1,  and  for  1^^   the  simplification  is  ar  «  ^.   This 
is  the  justification  of  equation  (4.17).   As  a  result 

oCC, 
I  ^ 


120 


And  for  Ir-,-,,    the  result  Is 


aT- 


I^^|-  <  O(a^) 


J 

ax 


^^{i)di   <  y^   M^f^^(r)  , 


which  Implies,  together  with  the  result  for  lypj 


A 


1^1  <  2  /a  M^fi^(r) 


(VI. 7) 


From  the  expansion  (4.24)  for  Q.{r,a)    It  Is  seen  that 


A 


cD^(r)  -  fi(r,a)|  <  /a  M^n^(r)  . 


(VI. 8) 


Summing  up  (VI.5),  (VI. 6)  and  (VI. 7)  gives 


|E(r,a)|  <  4/a  M-j^Q^(r)  ,     for   r  <  t, 


and  hence 


E(r,a)j  <  )/aM  ^'(r,a) 


for   r  <  T 


(VI. 9) 


Combination  of  (VI. 5)  and  (VI. 9)  Implies 


|E(r,a)|  <  /a  M  f2(r,a)      for   1  <  r  <  oo  , 


which  is  what  is  to  be  proved 
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